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ABSTRACT 

Lumped-parameter  and  finite-element-models  are  devel- 
oped for  the  dynamic  behavior  of  plane  straight  sided 
frames.   The  models  do  not  include  axial  and  shear  defor- 
mation, but  the  equations  of  motion  developed  using  the 
models  allow  for  time  varying  external  loading.   The  per- 
formance of  these  models  is  evaluated  by  a  comparison  with 
a  standard  transfer  matrix  method  for  the  special  case  of 
free  undamped  vibration.   The  finite-element-model  proves 
to  be  much  the  better  model.   For  the  first  five  modes, 
the  finite-element-model  with  27  degrees  of  freedom  differs 
by  no  more  than  0.4%  from  the  "exact  values"  given  by  the 
transfer  matrix  method.   The  lumped-parameter-model  with 
39  degrees  of  freedom  gives  errors  roughly  ten  times  as 
great. 
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NOMENCLAUTRE 


A.  =  Cross  sectional  area  of  member  i. 

B  =  Spectral  matrix  from  diagonalization  of  J. 

C  =  Matrix  relating  d  to  c. 

c  =  Coefficient  vector. 

c.  =  Element  of  coefficient  vector 

D  =  Modal  matrix  from  diagonalization  of  J. 

d  =  Deflection  and  slope  vector. 

E.  =  Modulus  of  elasticity  of  member  i. 

EI  =  Flexural  rigidity 

f_  =  External  force  vector. 

f  =  Vector  of  x  components  of  external  force. 
-— x 

f  .  =  x  component  of  external  force  at  node  j . 

f  =  Vector  of  y  components  of  external  force. 

— y 

f  .  =  y  component  of  external  force  at  node  j. 

t  d2w 

g_  =  Row  vector  in  powers  of  r  for  — »•  . 

dr^ 

h  =  Row  vector  in  powers  of  r  for  4>, 

I_  =  Identity  matrix. 

I.  =  Second  moment  of  area  A.  with  respect  to  a  centroidal 
axis  perpendicular  to  plane  of  frame. 

J  =  Matrix  sum  of  inertia  terms. 

J  =  Inertia  matrix. 
— o 

J.  =  Mass  moment  of  inertia  of  segment  j  about  an  axis 

-1  normal  to  the  plane  of  the  structure  and  thru  the 
segment  centroid. 

J. .  =  Rotational  inertia  between  j  and  j+1. 
33 

K  =  Stiffness  matrix. 


k.  =  Stiffness  of  segment  j. 

L  =  Incidence  matrix. 

s  . 

i.  =  Length  of  segment  of  member  i,  =  —  . 
1  n . 

1 

M  =  Mass  matrix. 

M  =  Rotational  mass  matrix. 
— r 

M  =  Translational  mass  matrix. 

m.  =  Segment  mass. 

J2»  =  Moment  vector. 

"lj  =  Moment  applied  between  node  j  and  j-1. 

N  =  Total  number  of  subdivisions. 

n.  =  Number  of  subdivisions  in  member  i. 

P  =  Matrix  used  to  apply  end  conditions  in  lumped- 
parameter-model  . 

£  =  Row  matrix  in  powers  of  r  for  w. 

Q  =  Moment  arm  matrix  for  external  force  vector. 

£  =  Generalized  force  vector. 

R  =  Corner  transformation  matrix. 

r  =  Axial  coordinate. 


2r.tT.  ^_  2 


S  =  Matrix  product  of  B  2D"K  DB 

s.  =  Length  of  member  i. 

T  =  Kinetic  energy. 

T  =  Rotational  kinetic  energy. 

T  =  Translational  kinetic  energy. 

U  =  Potential  energy. 

u  =  Vector  of  x  components  of  displacement 

u  =  Axial  displacement. 
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u.  =  x  displacement  of  node  j. 

v  =  Vector  of  y  components  of  displacement, 

v.  =  y  displacement  of  node  j. 

W.  =  Specific  weight  of  member  i. 

w  =  Transverse  deflection. 

X  =  Matrix  of  the  x  components  of  distance. 

X .  =  x-coordinate  of  the  centroid  of  segment  j . 

x. .  =  x  component  of  distance  from  node  j  to  node  i+1. 

Y  =  Matrix  of  the  y  components  of  distance. 

Y.  =  y-coordinate  of  the  centroid  of  segment  j. 

y. .  =  y  component  of  distance  from  node  j  to  node  i+1. 

*      t 

z  =  product  of  B   and  D 

(x,y,z)  =  Cartesian  coordinate  system. 

a  =  Deviation  angle  at  a  corner  of  the  frame. 

3  =  cos  a. 

Y  =  Angle  between  axial  direction  in  local  coordinate 

system  and  x  direction  in  Cartesian  system. 

n  =  sin  a. 

£  =  Relative  rotation  vector. 

0.  =  Relative  rotation  of  hinge  j. 

y  =  Mass  per  unit  length. 

<j>  =  Slope. 

w  =  Natural  circular  frequency. 

Subscripts 

a  =  assembled  matrices. 

i  =  Property  of  member  i. 

j  =  Property  at  node  j . 
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ij    =  Property  from  node  i  to  node  j  in  finite-element- 
model;  from  node  j  to  node  i+1  in  lumped-parameter- 
model. 

u    =  Upper  rows  of  a  matrix  as  determined  by  end  con- 
ditions. 

i  =   Remaining  rows  of  matrix. 

Superscripts 

t    =  Transpose  of  vectors  and  matrices. 

=  Matrix  after  application  of  end  conditions. 

*    =  Vector  resulting  from  corner  transformation. 

1     =  In  lumped-parameter-model  denotes  a  once  reduced 

matrix;  in  finite-element-model  denotes  the  proper- 
ties of  the  element  to  the  right  of  the  corner. 
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INTRODUCTION 

When  designing  a  structure,  the  engineer  must  be  aware 
of  the  many  dynamic  disturbances  that  may  act  upon  this 
structure.   He  must  also  be  aware  of  the  response  that  the 
structure  will  have  to  these  disturbances.   The  calculation 
of  the  response  is  usually  a  very  difficult  undertaking. 
To  assist  in  this  calculation,  the  structure  is  generally 
characterized  by  a  simplified  model,  and  the  analysis 
carried  out  using  this  model. 

In  this  paper,  two  models  for  analyzing  plane  frame 
structures  with  "chain"  topology*  are  developed.   The 
first  model  lumps  both  mass  and  flexibility  and  the  second 
is  a  finite-element- model.   To  test  these  models,  analyses 
to  find  the  in-plane  circular  frequencies  and  translational 
displacements  are  carried  out  on  a  typical  gable  bent.   The 
results  of  these  analyses  are  compared  to  frequencies  and 
displacements  obtained  through  the  use  of  the  so-called 
transfer  matrix  method. 

The  equations  of  motion  as  developed  for  the  models 
provide  for  time  varying  external  loads,  but  do  not  include 
damping.   The  analysis  of  the  gable  bent  is  carried  out  for 
the  special  case  of  free  vibration.   This  simplification 


*Plane  framed  structures  having  2  ends  and  a  single 
path  from  end  to  end. 

**Superscripts  refer  to  the  bibliography  which  appears 
on  page  48. 
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is  required  so  that  the  principal  modes  obtained  through 
the  use  of  the  models  may  be  compared  with  values  from  the 
transfer  matrix  method. 
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DEVELOPMENT  OF  THE  LUMPED-PARAMETER-MODEL 

In  order  to  find  in-plane  frequencies  and  displace- 
ments of  a  structure  consisting  of  straight  sides  with  uni- 
form cross  section,  a  model  lumping  both  mass  and  flexibil- 
ity is  developed.   This  model  neglects  axial  and  shear  de- 
formation, but  can  include  rotatory  inertia. 

The  model  is  formed  by  subdividing  each  straight  mem- 
ber into  segments.   These  segments  are  considered  to  have 
square  ends.   The  length  along  the  centerline  of  the 
straight  member  and  the  sum  of  the  lengths  of  the  segments 
replacing  it  are  the  same,  even  though  this  would  cause  some 
overlapping  at  the  corners.   The  mass  and  flexibility  of 
each  segment  is  lumped  at  the  segment  centroid.   The  struc- 
ture is  then  replaced  by  a  chain  of  straight  rigid  bars , 
pin  connected  at  the  locations  of  the  segment  centroids. 
The  mass  of  the  segment  remains  concentrated  at  the  hinge, 
and  the  flexibility  of  the  segment  is  replaced  by  a  tor- 
sional spring.   The  equations  of  motion  are  developed  and 
the  model  analyzed  to  obtain  its  natural  frequencies  and 
the  displacements  of  the  hinges. 

Rigid  body  translation  is  prevented  by  taking  the  left 
end  as  fixed  or  pinned.   The  right  end  is  allowed  to  assume 
any  of  the  common  end  conditions:  fixed,  pinned,  simply 
supported,  or  free. 
Development  of  the  Equations  of  Motion 

To  illustrate  the  method,  a  structure  containing  two 
straight  members  intersecting  at  an  acute  angle  is  used. 
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(a)  Two  Member  Frame 


A 


t=z=*~ 


(b)  Segmented  Frame 


Y7  _ 


t'^ 


T-^f 


x7 


(c)  Lumped  Model  of 
Frame 


Figure  1 

Both  ends  of  the  structure  are  taken  as  fixed  (Figure  la) . 
For  each  straight  member  i  (1=1,2),  let 
W.  =  specific  weight,  lb/cu.in. 

A.  =  cross  sectional  area,  sq.in. 

I .  =  second  moment  of  the  cross  sectional  area  with 
respect  to  the  z  axis,  in.14 

E.  =  modulus  of  elasticity,  lb/sq.in. 


s.  =  length  of  member,  in. 


n.  =  number  of  segments;  n,=2,  n2=3. 
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Each  section  of  the  structure  is  divided  into  n.  seg- 
ments (Figure  lb).   The  n.  need  not  be  the  same  for  the  dif- 
ferent members,  but  in  each  member  the  segments  are  of  uniform 
length.   The  parameters  of  each  segment  are  then  lumped  at  its 
centroid.   These  parameters  are 

I .    =  segment  length  =  s./n.  inches. 

2 

m.  =  segment  mass  =  (AWil)  ./386  lb-sec.  /in. 

J.  =  mass  moment  of  inertia  about  an  axis  normal 

to  the  plane  of  the  structure  and  through  the 
segment  centroid, 

=  m. (£2/12  +  I. /A.)  in-lb-sec? 

k.  =  segment  stiffness  =  (EI). /I.    lb-in. 

The  structure  is  then  replaced  by  the  rigid  bars  and  tor- 
sional springs  (Figure  lc) .   The  springs  are  unstrained  in  the 
datum  position.   The  Cartesian  coordinates  of  the  end  points 
and  hinges  of  the  structure  in  the  datum  position  are  taken 
as  X.  and  Y.  where  j  =  l  ,2  ,  3  ,4  , 5  ,6  ,  7 .   Thus  the  displacements 
of  the  hinges  (u.,  v.)  in  the  positive  x  and  y  directions  for 
a  relative  rotation,  0.,  at  hinge  j  are  given  by 


where 


u  =  -Y  0_ 
v  =   X  6 


u  =  col(u2,  u3,  u4,  u5,  u6,  u?) 
v  =  col(v2,  v3,  v4,  v5,  v6,  v?) 


e  =  coi(e1,  e2,  e3,  e4,  85,  e6) 


(1) 


17 


Y  = 


*51  Y52  Y53  Y54  Y55  ° 


Pel  Y62  Y63  Y64  Y65  Y66 


X  = 


xi;l  0   0   0   0   0 


x21  x22  0    0 


X31  X32  X33  °    °    ° 


X41  X42  x43  x44  ° 


X51  X52  X53  X54  "55 


! X61  X62  x63  X64  X65  X6 


0 


x55  °   1 


J 


y.  .  =  Y.  "  -  Y. 

Ji]     l+l     j 

x.  .  =  X.^,  -  X. 
13     l+l     j 

The  following  notation  has  been  used  in  this  equation: 

1.  Column  vectors  are  denoted  by  lower  case  symbols. 

2.  Rectangular  matrices  are  denoted  by  capitalized  symbols. 

3.  All  vectors  and  matrices  are  underlined. 

This  notation  will  be  used  in  all  other  such  equations. 

Associated  with  each  bar  is  a  rotational  inertia  made  up 
of  half  the  mass  moments  of  inertia  of  the  preceding  and  fol- 
lowing segments.   This  rotational  inertia  is  represented  by 
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J.  .  =  i(J  .  +  J.  .,)  . 

33      2     3        :+i 

Using  equation  (1),  the  kinetic  energy,  T,  and  the  poten- 
tial energy,  U,  of  the  structure  are  formed. 


where 


l»tt   •  1 • t  t   •    1 • t  t 

T  =  i0  X  MX0  +  ie  Y  MY6  +  i-6  L  J  L6 

2—  —  2—  —  '      2—   —  — o — 

U  =  ho^KQ 


M  =  diag(m2,  nu,  m .  ,  m- ,  mg  ,  0) 

^o=  diag(JU/  J22'  J33'  J44'  J55 '  J66) 
K  =  diag(k1,  k2 ,  k3,  k4 ,  kg  ,  kg) 


L  = 


1 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

1 

1 

1 

0 

0 

0 

1 

1 

1 

1 

0 

0 

1 

1 

1 

1 

1 

0 

1 

1 

1 

1 

1 

1 

(2) 


Derivatives  with  respect  to  time  are  denoted  by  dots  and  the 
transpose  of  a  matrix  by  a  superscript  t. 
The  Lagrangian  equations  of  motion, 


may  be  written 


where 


dt 


3T\_  8T    8U  , 


39 


86    39 


J9  +  K9  =  q 


J  =  XfcMX  +  YfcMY  +  LfcJ  L 


q  =  generalized  force  vector. 
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(3) 


Let  the  x  components  of  the  external  forces  applied  at  hinges 

2  to  7  be  represented  by  the  vector 

f  =  col(f  J.    f'  -j ,  f  „  .  f  c ,  f  * ,    f  _,)  • 
—x         x2 '   x3'   x4 '   x5 '   x6    x7 

In  a  like  manner  define  a  vector  f  .   For  the  counterclockwise 

-y 

external  moments  Yt[ .  applied  between  hinge  i-1  and  i,  the  vec- 
tor 

h^=  col(*i2,  K3/  *l4'  ^5^  >n6/  hX7) 

is  defined.   The  generalized  force  vector  may  be  expressed  in 
terms  of  these  vectors  as 

a.=  (iV"  Y%  +  ^%  •  (4) 

It  should  be  noted  that  the  matrices  J  and  K  are  both 
symmetric.   This  condition  results  when  any  symmetric  matrix 
is  postmultiplied  by  a  matrix  and  is  premultiplied  by  its 
transpose.   This  symmetry  property  is  highly  desirable  when 
dealing  with  matrix  eigenvalue  problems  and  should  be  preserved 
if  possible  when  applying  boundary  conditions. 
Application  of  Boundary  Conditions 

Up  to  this  point  the  development  has  been  general  in  na- 
ture.  We  now  restrict  further  development  to  free  vibration. 
With  this  restriction,  the  only  external  forces  applied  are 
those  required  to  satisfy  end  constraints.   If  the  model  has 
both  ends  fixed,  the  end  conditions  and  corresponding  exter- 
nal forces  are 

e1  =  o,  *n2  ft   o 

v?  =  0,  fy7  t    0 
u7  =  0,  fx7  *   0 

e1  +  e2  +  e3  +  e4  +  e5  +  e6  =  o,  *n7  f   o  . 
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Equation  (4)  can  therefore  be  reduced  to 

a   =  Q  1 


(5) 


where 


f  =  col(r*i2,  fy?,  -  fx?,  K-j) 

a  x. 


:61   y61 


0   x 


0   x 


62   y62 


Q  = 


63   *63 


0   x 


64   y64 


0   x 


65   y65 


0   x 


66   y66 


The  columns  of  Q  are  the  first  and  last  columns  of  L  ,  and 

t      t 
the  last  columns  of  X   and  Y  . 

The  end  conditions  may  be  represented  by 


ere  =  0  . 


(6) 


This  equation  is  derived  from  the  following  set  of  equations. 
8.  =  0 

u7  =  x61e1+x62e2+x63e3+x64e4+x6505+x6606  =  o 

V7  =  y6iei+y62e2+Y63e3+y64e4+y65e5+Y66e6  =  ° 

e,  +  e0  +  e0  +  e„  +  ec  +  er   =  o 

1     2     i     4     5     6 
A  process  must  be  devised  to  incorporate  the  boundary 
conditions.   This  process  will  reduce  the  number  of  degrees  of 
freedom  by  an  amount  equal  to  the  number  of  constraint  condi- 
tions, in  this  case  four.   The  result  will  be  a  pair  of  homo- 
geneous equations  with  two  remaining  6.'s. 
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A  normal  Gaussian  elimination  technique  for  the  solution 
of  simultaneous  equations  is  used  to  solve  for  f  _,  f  _,  >r\~ , 
and  yy\^.      Through  the  use  of  the  end  conditions,  the  number  of 
unknown  6  .  is  next  reduced  to  two.   The  four  non-homogeneous 
equations  containing  the  forces  and  moments  can  be  deleted, 
leaving  only  the  two  homogeneous  equations.   A  similar  elim- 
ination could  be  used  for  any  other  set  of  end  conditions. 

This  process  can  be  represented  by  a  set  of  matrix  opera- 
tions.  Equation  (3)  is  written  in  partitioned  form  as 


(7) 


where  the  subscript  u  represents  the  rows  which  will  contain 
the  non-homogeneous  equations,  and  I   the  remaining  rows.   This 
may  require  a  row  interchange  to  insure  a  non-singular  Q  . 
In  the  present  case  Q   is  4x4.   Gaussian  elimination  is  now 
performed  to  solve  for  the  forces  and  moments.   In  this  proc- 
ess the  elements  of  Qg   will  be  reduced  to  zero.   The  lower 
portion  of  equation  (7)  may  be  rewritten 


J 

— u 

— u 

"au~ 

A 

£  t 

.-* 

£  ■ 

ft* 

I'^i   +  k  J  e  £  =  o 


(8) 


where  the  prime  denotes  reduced  matrices. 

The  end  conditions,  equation  (6),  may  be  rewritten 


t   t 


0 
— u 


-I 


=  0 


and  from  this  we  get 


Qt6   +  Q^60  =  0 
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and  therefore 


6   =  -(Qt)~lQi:Q() 
— u     — u   —  I— I 


Thus  _6   may  be  eliminated  from  _6  by  the  transformation 


~<^t 


h   =  *i* 


(9) 


where  I_  is  the  identity  matrix,  in  this  case  of  order  2x2. 
It  is  easily  shown  that 


PfcQ  =  0 


and  that  in  equation  (8) 


and 


ii  =  *-tj- 


£i  =  %-&  - 


The  Q_   may  then  be  replaced  in  equation  (8)  by  P9.  .   After  dif- 
ferentiating equation  (9)  twice,  the  8_  of  equation  (8)  may 
also  be  replaced,  and  the  following  expression  obtained: 


where 


llh   +  K»e  =  o 


j"  =  pfcjp 


K£  =  PUKP 


(10) 


The  boundary  conditions  can  therefore  be  applied  and  a  set  of 
homogeneous  equations  obtained  by  this  process.   This  trans- 
formation does  not  destroy  the  symmetry  of  the  matrices. 
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Eigenvalue  Solution  of  the  Equations  of  Motion 

Assuming  that  the  vibration  will  be  sinusoidal,  equation 
(10)  may  be  written 

Elh  -  -2llh  (ID 

where  a)  is  the  natural  circular  frequency.   This  equation  is 
in  a  form  suitable  for  eigenvalue  solution. 

The  solution  of  equation  (11)  is  now  undertaken.   This 

2 

solution  will  involve  eigenvalues  and  eigenvectors.    A  sym- 
metric matrix  may  be  expressed  in  terms  of  its  modal  matrix 
and  its  spectral  matrix  where  the  spectral  matrix  is  a  diag- 
onal matrix  of  eigenvalues  and  the  columns  of  the  modal  matrix 
are  eigenvectors.   This  property  is  used  several  times  in  the 
solution.   A  further  important  property  of  the  modal  matrix 
is  that  the  eigenvectors  are  normalized  with  respect  to  the 
identity  matrix.   Thus  the  transpose  of  the  modal  matrix  is 
also  its  inverse. 

If  J"  is  resolved  into  its  (diagonal)  spectral  matrix 
B  and  modal  matrix  D  according  to  the  relation 


J"  =  DBD' 


equation  (11)  becomes 


K£_6£    =    oa2DBDt_6£    . 

Since   DD      =    I_,    the    above   equation  may   be   rewritten 

KJDD*^    =    w2DBDt_9£    . 

Since  B  is  a  diagonal  matrix  of  positive  elements,  we  define 

i 
the  principal  square  foot  as  B2  where  the  elements  of  this  ma- 
trix are  the  positive  square  roots  of  the  elements  of  B.   Not- 

_2_  l 
ing  that  B  232  =  I,  the  following  form  results, 
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K^DB  2B_2DU6_   =  ai  DB2B2D  8,  . 


I  t 
Letting  z_  =  B_2p_  fl.  , 


i      ->   i 


2 
Solving  for  u>  z, 


KJDB;  2z  =  w2DB2z  . 


-it    --J     2 

B  2D  K"DB  2z  =  to  Z 


-i  t     -i 
Letting  S  =  B  2D  K"PB  2,  the  equation  may  be  rewritten 

2 

SZ  =  CO  £  . 

2 

By  diagonalization  of  S_,  the  values  of  w   and  z   may  be  ob- 
tained.  From  these  values,  the  natural  frequencies  are  ob- 
tained directly.   The  values  of  the  rotations  of  the  nodes 

t 
are  obtained  from  the  equation  6.  =  DB2z_.   The  rotations 

at  all  the  hinges  of  the  structure  can  be  recovered  by  use 
of  equation  (9) .   The  translation  displacements  of  the  hinges 
are  easily  obtained  with  the  use  of  equations  (1) . 
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DEVELOPMENT  OF  FINITE-ELEMENT-MODEL 
General 

A  finite- element- model  is  developed  to  describe  the 
in-plane  motion  of  a  frame  with  straight  members  of  uniform 
cross  section.   The  model  neglects  axial  and  shear  deforma- 
tion and  can  include  rotatory  inertia  if  desired. 

The  use  of  a  finite  element  approach  requires  the  divi- 
sion of  the  structure  into  elements «  and  the  determination 
of  element  mass  and  stiffness  matrices.   The  determination 
of  the  elemental  mass  and  stiffness  matrices  will  be  taken 
up  first. 

Consider  a  general  element  of  length  1,  in  a  local 
coordinate  system.   In  this  system,  r  is  the  axial  coor- 
dinate, w  is  the  transverse  deflection,  and  $  is  the  slope. 
The  element  has  no  flexibility  in  the  r  direction.   If  we 
assume  that  w  may  be  approximated  by  a  third  degree  polyno- 
mial in  r,  then 

2       3 

w  =  c  +  c,r  +  c~r  +  c-.r 
o     1      2       3 


and 


dw        ~       -    2 
d>  =  3—  =  c,  +  2c~r  +  3c0x 
Y    dr     1      2       3 


where  c  ,  c, ,  c2 ,  c~  are  unknown  coefficients.   These  equa- 
tions may  be  expressed  as 

t 
w  =  p_  c 


<J>  =  hfcc 


where 


t  2    3 

p_  =  row(l,  r,  r  ,  r  ) 

t  2 

h   =  row(0,  1,  2r,  3r  ) 


£  =  col(c  ,  c,  ,  c2 ,  c^)  . 


26 


If  w  and  <j>  are  evaluated  at  the  ends  of  the  element,  the  re- 
sults may  be  expressed  as 


where 


d  =  C"1c 


d  =  col(wQ,  <J>0,  w£,  <bz) 


10    0    0 
0    10    0 


-1 


1    I         I2      I3 


0    1  21       3£: 


and  the  subscripts  denote  evaluation  at  r=0 ,  and  r=l.      By 


-1 


inverting  C   we  may  solve  for  c, 


where 


c  =  Cd 


0    0    0 


(12) 


\  I2    I 


0    0 
3   -1 


_2   1   -2 
£3   £2  £3 


Bending  Energy 

The  bending  energy,  U,  of  the  element  is  given  by 


°-n 


k/4-51  dr 


(13) 


'0   \dr 
where  EI  is  the  flexural  rigidity  for  in-plane  bending. 

Now 


d  w    „      , 
— =-  =  2c0  +  6c0r  =   t 
2      3    £  c 


where 


dr 


g   =  row(0,  0,  2,  6r)  . 
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Since  EI  is  a  scalar  constant,  the  equation  may  be  rewritten 


1  (     t   t 

U  =  i-EI  \  cugg_  cdr  . 


Since  c  and  c   are  constants,  they  may  be  removed  from  the 
integral.   Substituting  equation  (12)  for  c  and  c  we  may 
write 

U  =  ^K'd  (14) 

where 


K'  =  bending  stiffness  =  EICt  ( Cg_£tdr)  C  .     (15) 

It  can  be  shown  by  carrying  out  the  required  integration  and 
matrix  multiplications  that 


K'  =(EI/l3) 


12 

61 

-12 

61 

61 

4£2 

-  eiz 

2£2 

-12 

-61 

12 

-61 

en 

2£2 

-  61 

4£2 

(16) 


Kinetic  Energy 

The  kinetic  energy,  T,  of  the  element  is  given  by 


■  H"2dr  +  5^Sn*2dr       (l7) 


where 


y  =  mass  per  unit  length 

A  =  cross  sectional  area 

I  =  second  moment  of  area  A  about  an  axis  through 
the  centerline  of  the  element  and  perpendicular 
to  both  the  w  and  r  axes • 

The  first  term  is  the  translation  energy,  T  ,  and  the  second 

t 

term  is  the  rotational  energy,  T  .   Thus 


i  (l 


w  dr 
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and 


m     1  ul C   i2  . 


Now 


2    •   3 
w  =  c~  +  c,r  +  c-r  +  c~r 

t-     t  * 
=  £  c  =  £  Cd  . 

Again  noting  that  C  and  C   do  not  depend  on  r,  T  may  be  re- 


written 


where 


Using 


*t  =  H*s&s 


**t  =  ^-iS  E£-tdrj-  • 


4>  =  c,  +  2c~r  +  3c^r' 


t  •     t  ' 
=  he  =  h  Cd  , 


by  a  similar  process 


T   =  id^'d 
r   2 r— 


where 


Si-^c^^c, 


After  evaluation 


M.'  =  (y£/42  0) 


156     22£     54     -13£ 
22£     4lz         131        -    3V 
54     13£    156     -22£ 


-  13£   -  3£2  -  22£ 


4£ 


and 


M'  =(yI/30AX) 


36      31  -  36  35~ 

3£     4£2  -  31  -I2 

-  36  -  3£  36  -  31 

3£  -   £2  -   3£  4£2 


(18) 


(19) 


(20) 


(21) 


22) 


(23) 
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If  equations  (18)  and  (20)  are  recombined,  then 


where 


T  =  ^d^M'd  (24) 


M'  =  M'  +  M1  (25) 

—    — t   — r 


3 

and  M*  is  called  the  consistent  mass  matrix  for  bending. 

Axial  Displacement 

Up  to  this  point  there  has  been  no  mention  of  a  dis- 
placement in  the  axial  direction.   Since  the  element  has  no 
axial  flexibility,  any  displacement  in  the  axial  direction 
will  be  displacement  of  the  element  as  a  rigid  body.   This 
displacement  will  be  the  same  at  either  end  of  the  element. 
The  axial  displacement  is  represented  by  u,  where  u=u  =u„ . 
In  fact,  as  will  be  shown  later,  the  axial  displacement  of 
an  entire  straight  member  made  up  of  several  elements  may 
be  represented  by  a  single  symbol. 

The  axial  displacement  should  be  included  in  the  bending 
and  kinetic  energy  equations.   The  u  has  no  effect  on  the 
bending  energy  or  rotational  energy.   It  does  contribute  to 
the  kinetic  energy  of  translation.   This  contribution  is 

-i       <y 

yyilu    .      Equations    (14)    and    (24)    may   be   rewritten 

U  =  i  y^Kv  (26) 

T  =  \  xlix  (27) 


where 


v  =   col(u,   wQ,    cj>0,    w^,    cj^) 


K  = 


0! 

0 

0 

0 

0 

0! 

0! 

o: 

K' 

0' 

(29) 
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M  =  M.  +  M 

—   — t   — r 


(30) 


_o_ 

_0__0_ 

_0_ 

«t= 

o   ! 

o  : 
o  ; 

— t 

^o  ; 

0 

0      0 

0  ~~ 

o  ; 

M    = 
— r 

0    J 

o  : 
o  ; 

M' 
— r  ■ 

Corner  Transformation 

The  combination  of  the  element  kinetic  and  bending  en- 
ergy to  form  the  total  bending  energy  and  kinetic  energy  of 
the  structure  is  not  difficult,  except  at  a  corner.   Figure 
2a  shows  two  elements,  one  to  the  left  and  one  to  the  right 
of  the  corner.   A  subscripting  notation  has  been  added.   A 
single  subscript  is  used  to  denote  the  end  or  node  at  which 
the  quantity  is  evaluated.   A  double  subscript  is  used  to 
denote  quantities  which  are  the  same  at  each  end  of  an  el- 
ement, and  for  vectors  and  matrices  which  pertain  to  the 
element.   Thus  v,2  =  col(u,2,  w,  ,  $,,  w2 ,  <f>2)  /  and  so  on. 


3   rr' 


(a)  Elements  Adjacent  to 
Corner 


(b)  Coordinates  at 
Node  2 


Figure  2 
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The  displacements  at  node  2  are  with  reference  to  two 
different  local  coordinate  systems,  and  cannot  be  directly 
combined  (Figure  2b) .   The  prime  is  used  to  denote  the  dis- 
placements at  node  2  which  are  in  the  local  coordinate  system 
of  the  element  to  the  right  of  the  corner.   There  are  only 
two  independent  displacements  at  node  2.   It  is  thus  neces- 
sary to  express  the  four  displacements  at  this  node  in  terms 
of  only  two  displacements.   The  axial  displacements  of  the 
two  elements  are  chosen.   This  choice  is  desirable  since  the 
axial  displacement  of  all  parts  of  a  straight  member  is  the 
same.   Further,  when  end  conditions  are  applied,  the  axial 
displacements  at  the  ends  are  some  of  the  quantities  which 
are  known . 

Letting  3  =  cos  a  and  n  =  Sin  a,  from  figure  2b  we 

can  write 

w2,  =  3w2  -  nu12 

u2,3  =  nw2  +  3u12  . 

Using  the  above  equations,  w2  and  w2 ,  may  be  expressed  in 
terms  of  u,2  and  u.,,  as: 

w2,  =  -(l/n)u12+  (3/n)u2,3  (31) 

w2  =  -(B/n)u12  ■■■+  (l/n)u2,3         (32) 

This  transformation  may  be  put  into  a  convenient  matrix 
notation.   For  the  element  to  the  left  of  the  corner,  this 
transformation  is  written 

X12  =  5Xl2 
1Z  LZ  (33) 

-i2  -   <Z&V 
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where 


v*2  =  col(u12,  wx,  cf)1,  u2,3,  <p2) 


R  = 


1 

0 

0 

0 

0 

o 

1 

0 

0 

0 

0 

0 

1 

0 

0 

-I 

n 

0 

0 

1 

n 

0 

o 

0 

0 

0 

i 

The  kinetic  and  potential  energy  of  the  element  may  there- 
fore be  written 


T12  -   j(Xi2)VM12RvJ2 


U12  =   7(X12>VKi25v!2   ■ 


(34) 


For   the    element   to    the    right   of    the    corner,    the    transforma- 
tion  is 

v0  ,  -j    =   R   vi  ,  ,, 


where 


-2' 3         - 


v2,3    -     (v*,3)t(R-)t 


v*,3   =   col(u12,    u2,3,    <J>2,    w3,    <J>3) 


(35) 


R'    = 


1 

n 
o 
o 

0 


0 

0 

0 

0 

0 

0 

1 

0 

0 

1 
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The  kinetic  and  bending  energy  for  the  element  to  the  right 
of  the  corner  are  written 


T2'3  =  I^S^'^^R'v*^ 


u2'3  =  T(2S.3)t(£,)tTE2.3S,Z!-3  • 


(36) 


The  kinetic  and  bending  energies  in  equations  (34)  and  (36) 
can  now  be  combined. 
Example  of  Method 


(a)  Two  Member  Frame 


(b)  Frame  Divided 
Into  Elements 


Figure  3 

To  illustrate  the  method,  a  frame  consisting  of  two 
straight  members  intersecting  at  an  acute  angle,  a,  is  used 
(Figure  3a) .   Both  ends  of  the  structure  are  fixed. 

Each  member  is  divided  into  2  elements  (Figure  3b) . 
Each  of  the  two  members  need  not  have  been  divided  into  the 
same  number  of  elements.   The  nodes  at  the  ends  and  between 
elements  are  numbered  for  identification.   The  elemental  mass 
and  stiffness  matrices  are  determined  and  the  following  set 
of  equations  written: 

T12  =  fei2^12Xl2  '  U12  =  kl2*12Xl2         (37) 
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T2  3  =  ^2  3^2  3^2  3  '  U2  3  =  1X2  3^2  3X2  3         (38) 


l.'.t  „   _•_      TT      It 

-3< 

t 


T34  =  2X3A4X34  '  U34  =  2X34^34X34         (39) 


T45  "  1X45^45X45  '  U45  =  1X45^45X45   '      (40) 

The  corner  transformation  must  be  performed  at  node  3. 
Equations  (38)  and  (39)  become 

T23  "  I'ila3 ^23^23  '  U23  =  l<2&>  V^^ls   (41) 

T3'4=  |(i3.4»t(R')tll3'4£'i3'4  ' 

(42) 

The  elemental  kinetic  and  bending  energy  equations  are  com- 
bined and  are  written 


1*  t 
T   =  iv  M  v 
a   2— a— a— a 


U  =  ^vtK  v 
a   2— a— a— a 


(43) 


where 

y_a  =    col(u13,    w1,    <J>^,    yj2,    <J>2,    u35,    <j>3,    w^,    cj>4,   w5 ,    <|>5) 

M   =  assembled  mass  matrix 
—a 

K   =  assembled  stiffness  matrix 
—a 

and  the  subscript  a  denotes  assembled  quantities.   The  assem- 
bled matrices  are  formed  by  placing  the  elemental  matrices 

in  their  correct  positions  in  the  larger  framework  of  the 

4 
assembled  matrix,  and  summing  all  the  overlapping  terms. 
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The  Equations  of  Motion 

Using  the  Lagrangian  equation  of  motion, 

equation  (43)  may  be  written 

M  v   +  K  v   =  q  (44) 

—a— a   —a— a    ■* 

where  q_  is  the  generalized  force  vector.   The  elements  of  q_ 

are  the  generalized  forces  and  moments  corresponding  to  the 

elements  of  the  displacement  vector,  v  .   This  equation  is 

•■a 

general  in  that  no  restriction  has  been  placed  on  the  gen- 
eralized force  vector. 
End  Conditions 

We  now  restrict  further  development  to  free  vibration 
only.   Since  in  Figure  3  both  ends  are  fixed  and  there  is  no 
external  loading,  the  only  non-zero  forces  in  the  vector  q_ 
are  the  forces  required  to  satisfy  end  constraints. 

For  both  ends  fixed,  the  following  is  known  about  the 
structure : 

u13  =  w1  =  <p±   =    0 

u35  =  w5  =  d>5  =  0  . 

The  order  of  the  matrices  in  equation  (44)  is  now  reduced  by 

an  amount  equal  to  the  number  of  constraint  conditions,  in 

this  case  six.   The  rows  and  columns  corresponding  to  these 

conditions  are  deleted  and  the  equation  rewritten 

M"v"  +  K"v"  =  0  (45) 

—a— a   —a— a 
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where 

v^  =  col(w2,  <$>2,    cj)3,  w4,  (j)4) 


K"  =  a  5x5  reduced  matrix 
—a 


M"  =  a  5x5  reduced  matrix. 

- -  a 


Eigenvalue  solution 

Assuming  the  vibrational  motion  is  sinusoidal,  the 

equation  becomes 

K"v"  =  oo2M"v"  (46) 

—a— a     —a— a 

where  co  is  the  natural  circular  frequency.   This  equation  is 
now  in  a  form  suitable  for  eigenvalue  solution.   The  solution 
is  carried  out  in  the  same  manner  as  described  on  pages  24 
to  25  of  the  lumped  parameter  development.   This  solution 
will  yield  the  first  five  natural  frequencies  (since  matrices 
are  5x5) ,  the  transverse  displacements  at  the  2nd  and  4th 
node,  and  the  slopes  at  nodes  2,  3,  4.   The  displacements  are 
known  at  the  ends  from  the  end  conditions,  and  through  use  of 
equation  (31)  or  (32)  the  displacement  at  node  3  may  be  ob- 
tained.  These  displacements  are  all  referred  to  local  coor- 
dinate systems.   They  are  transformed  to  a  Cartesian  (x,y,z) 
system  through  the  use  of  the  following  equations: 

u.  =  -w .  sin  Y4  +  u.  cos  y. 
3  3  i'3       i 


where 


v.  =  w.  cos  Y-j'+  u-;  Sin  Yi 


Y-  =  angle  between  the  x  axis  and  axial 
direction  of  member  i 
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u.  =  displacement  in  x  direction 
v.  =  displacement  in  y  direction 
j  =  node  number,  j=l,2,3,4,5 
i  =  member  number,  i=l,2  . 
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PRINCIPAL  MODE  COMPARISON 


o 


-II 


8WF20 


30 


Figure  4 
Gable  Bent  Used  in  Analyses 

Figure  4  shows  the  structure  used  for  comparison  of  the 
models.   The  members  are  all  steel  wide-flanged  I-beams  of 
the  type  shown  on  the  figure.   Data  for  these  I-beams  is 
taken  from  reference  5. 
Transfer  Matrix  Method. 

The  transfer  matrix  method  is  widely  used  for  the  deter- 
mination of  the  principal  free  vibrational  modes  of  framed 
structures  and  will  be  used  as  the  standard  for  comparing 
the  finite-element  and  lumped-parameter-models.   The  devel- 
opment and  use  of  the  transfer  matrix  method  are  given  de- 
tailed treatment  in  reference  1.   A  computer  program  to  an- 
alyze similar  structures  was  developed  in  1964  by  Fink  in 
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his  master's  thesis  at  the  Naval  Postgraduate  School.    The 
"exact  values"  used  to  compare  the  models  are  generated 
through  the  use  of  a  transfer  matrix  computer  program  devel- 
oped by  the  author.   This  program  is  included  as  Appendix  D. 

The  first  five  undamped  natural  frequencies  and  corre- 
sponding mode  shapes  are  calculated  using  this  program  for 
four  end  conditions:   clamped-clamped,  clamped- simply  sup- 
ported, pinned-pinned,  and  pinned-simply  supported.   The 
frequencies  are  given  in  Table  I  and  the  mode  shapes  are 
listed  in  Appendix  A,  pages  50,  51,  52,  53.   In  Figure  5, 
the  mode  shapes  for  the  clamped-clamped  end  condition  are 
plotted.   The  displacements  are  greatly  exaggerated  so  that 
the  mode  shape  can  be  more  readily  visualized.   Only  the 
first  five  frequencies  are  tabulated,  since  after  this  the 
computer  round-off  error  begins  to  affect  the  mode  shape. 
The  frequencies  have  become  large  enough  that,  even  though 
the  convergence  routine  has  estimated  them  to  nine  decimal 
places,  the  frequency  determinant  is  not  zero.   Thus  when 
the  state  vector  of  the  left  end  is  passed  through  the  trans- 
fer matrices,  the  accumulation  of  round-off  error  is  evident 
at  the  right  end.   The  beginnings  of  this  error  can  be  seen 
in  the  computer  output  for  the  clamped-clamped  case  (page  50) . 
In  the  fifth  mode,  the  x  displacement  is  0.0014  at  the  right 
end,  but  should  be  zero  since  the  end  is  fixed.   This  sus- 
ceptibility of  the  transfer  matrix  method  to  round-off  error 
is  a  major  reason  why  other  computerized  models  are  sought. 

The  structure  was  turned  around  and  the  frequencies 
and  mode  shapes  calculated  for  the  pinned-pinned  and 
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clamped-clamped  end  conditions.   The  frequencies  are  un- 
changed, but  the  fifth  mode  shape  differs  slightly  in  the 
fourth  decimal  place.   This  difference  is  due  to  the  accu- 
mulated round-off  error.   The  outputs  for  these  runs  are 
listed  in  Appendix  D  on  page  112. 
Lumped- Parameter- Model 

An  analysis  of  the  same  gable  bent  using  the  lumped- 
parameter-model  is  next  carried  out.   Since  both  ends  clamped 
is  the  most  severely  constrained  condition,  this  end  condi- 
tion will  cause  the  structure  to  have  the  fewest  degrees  of 
freedom.   Using  this  end  condition,  it  is  expected  that  cal- 
culation of  the  natural  frequencies  and  displacements  will 
require  the  greatest  number  of  elements.   Therefore,  by 
successive  trials,  sufficient  subdividing  is  carried  out 
such  that  the  error  in  the  primary  frequency  is  less  than  1% 
and  the  error  in  the  5th  frequency  is  less  than  3%.   To  meet 
these  criteria  requires  38  subdivisions.   Using  the  same  num- 
ber of  subdivisions,  the  rest  of  the  end  conditions  are  run 
and  the  results  given  in  Table  I.   The  pinned-pinned  end  con- 
dition is  seen  to  lead  to  a  slightly  greater  per  cent  dif- 
ference than  the  clamped-clamped  condition. 

A  comparison  of  the  mode  shapes  for  this  model  with 
those  of  the  standard  is  made  by  comparing  the  outputs  on 
pages  50/  51/  52  > and  53  in  Appendix  A.   The  mode  shapes  agree 
to  two  decimal  places  for  the  first  four  modes,  and  to  only 
one  decimal  place  in  the  fifth  mode. 

The  error  due  to  the  lumping  of  the  parameters  decreases 
slowly  as  the  number  of  subdivisions  increases.   By  comparing 

41 


Plot  of  Mode  Shapes  for  a  Gable 
Pent  With  Both  Ends  Fixed. 


1st  mode 


2nd  mode 


\w 


3rd  mode 


4th  mode 


Figure  5 
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page  76  to  the  values  in  Table  I,  the  small  improvement  in 
modal  frequencies  resulting  from  increasing  the  number  of  sub- 
divisions can  be  seen.   There  is  a  similar  improvement  in  mode 
shapes.   If  an  accurate  estimate  of  frequency  and  displace- 
ment were  desired,  a  very  high  degree  of  discretization  would 
be  needed.   The  computer  time  necessary  to  solve  for  these 
values  would  be  large.   This  is  due  to  the  eigenvalue  routine 
where  the  time  for  solution  varies  as  the  cube  of  the  number 
of  subdivisions. 

The  frame  was  turned  around  and  values  obtained  for  the 
pinned-pinned  and  clamped-clamped  end  conditions.   These  fre- 
quencies and  displacements  (page  75)  are  identical  to  those 
calculated  in  Table  I.   This  duplication  of  results  adds  a 
measure  of  confidence  to  the  integrity  of  the  solution  using 
this  model. 
Finite  Element  Model 

The  frame  is  again  analyzed,  this  time  using  the  finite- 
element-model.   The  analysis  is  started  using  both  ends 
clamped.   To  meet  the  criteria  of  less  than  1%  difference 
in  the  primary  frequency  and  less  than  3%  difference  in  the 
5th  frequency  requires  12  elements.   The  frequencies  for  all 
end  conditions  using  12  elements  are  listed  in  Table  I. 

The  mode  shapes  obtained  by  this  method  are  in  excellent 
agreement  with  those  of  the  standard.   A  comparison  of  the 
mode  shapes  from  this  method  to  those  of  the  standard  (Appen- 
dix A,  pages  50,  51,  52,  and  53)  shows  agreement  to  three 
decimal  places. 
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TABLE  I 
Summary  of  Modal  Frequencies 

Clamped-Clamped  End  Condition 


Mode 
Number 

Transfer  Matrix 

Lumped-Parameter 
Sections : 6-13-13-6 

Finite-Element 
Sections :2-4-4-2 

Frequency 

Frequency 

%  dif. 

Frequency 

%  dif. 

1 
2 
3 
4 
5 

25.754 

34.680 

79.630 

150.284 

279.474 

i 

25.921 

34.929 

79.366 

148.262 

271.688 

+  0.65 
+0.72 
-0.33 
-1.34 
-2.78 

25.752 

34.679 

79.647 

150.444 

280.597 

-0.008 
-0.003 
+0.020 
+0.110 
+0.410 

Clamped-Simple  End  Condition 


Mode 
Number 


Transfer  Matrix 


Frequency 


Lumped-Parameter 
Sections : 6-13-13-6 


Frequency 


dif 


Finite-Element 
Sections:2-4-4-2 


Frequency  j%  dif. 


14.764 

34.236 

75.065 

123.283 

152.752 


14 

34 

75 

123 

152 


765 
238 
052 
234 
589 


14.769 

34.482 

74.760 

123.869 

150.914 


+  0.02 
+  0.71 
-0.39 
+0.51 
-1.11 


-0.008 
-0.006 
+0.018 
1+0.040 
+0.100 


Pinned-Pinned  End  Condition 


yiode 
Sf  umber 

Transfer  Matrix 

Lumped-Parameter 
Sections : 6-13-13-6 

Finite-Element 
Sections:2-4-4-2 

Frequency 

Frequency 

%  dif. 

Frequency ;%  dif. 

1 
2 
3 
4 
5 

15.257 

25.365 

75.652 

149.662 

277.069 

15.285 

25.365 

75.162 

147.574 

269.188 

+0.184 

0.00 
-0.65 
-1.39 
-2.84 

15.256  i-0.008 

25.364   -0.006 

75.667  j+0.020 

149.820  1+0.110 

278.160  '+0.390 

Pinned-Simple  End  Condition 


Mode 
Number 

Transfer  Matrix 

Lumped-Parameter 
Sections : 6-13-13-6 

Finite-Element 
Sections :2-4-4-2 

Frequency 

Frequency  1  %  dif. 

Frequency 

%  dif. 

1 
2 
3 
4 
5 

9.742 

24.419 

72.489 

122.907 

152.252 

9.747  j  +0.051 

24.419     0.00 

72.041    -0.62 

123.511  j  +0.49 

150.555    -1.11 

9.741 

24.417 

72.502 

122.956 

152.412 

-0.008 
-0.006 
+0.018 
+0.040 
+0.110 

44 


As  the  number  of  elements  used  increases,  the  differences 
between  the  calculated  values  and  the  standard  show  a  marked 
decrease.   The  increasing  of  the  subdivisions  from  8  to  12 
decreases  the  per  cent  difference  by  a  factor  of  10.   The 
results  for  the  8  element  structure  are  on  page  9  6  of  Appen- 
dix C.   This  marked  decrease  in  the  per  cent  difference  il- 
lustrates the  fact  that  this  model  will  give  good  results 
with  little  subdividing. 

It  is  known  that  the  finite-element-method  is  equivalent 

to  the  Rayleigh-Ritz  procedure  for  finding  natural  frequen- 

3 
cies.    Accordingly  the  estimate  for  the  frequencies  should 

represent  an  upper  bound.   The  values  in  Table  I  do  not 
agree  with  this.   In  all  cases,  the  first  and  second  frequen- 
cies from  the  finite-element-model  are  slightly  lower  than 
those  of  the  transfer  matrix  method.   An  attempt  to  find  a 
reason  for  this  discrepancy  in  the  programming  of  the  models 
was  unsuccessful. 

The  pinned-pinned  and  clamped-clamped  end  conditions 
were  rerun  with  the  structure  reversed.   The  frequencies  and 
mode  shapes  were  identical  to  those  previously  calculated 
(Appendix  C,  page  95). 

Using  the  finite-element-model  with  12  elements,  fre- 
quencies and  mode  shapes  were  calculated  with  rotatory  iner- 
tia included.   The  results  are  given  in  Table  II.   These 
results  show  that  rotatory  inertia  has  little  effect  on  the 
frequencies.   The  maximum  difference  between  frequencies 
with  and  without  rotatory  inertia  is  0.3%. 
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TABLE  II 
Effect  of  Rotatory  Inertia  on  Frequency 


!End 

— 

1 

Frequence 

{    (rad/sec) 

Condition 

1 

2 

3 

r 

4 

5 

1 
Clamped 

'Clamped 

W/0  ?l 
With 
%  dif. 

25. 

25. 

i  o. 

i 

752 
748 
015 

34.679 

34.675 

0.012 

79.647 

79.585 

0.078 

150. 

150. 

0. 

444 
207 
158 

280. 

279. 

0. 

597 
725 
312 

Clamped 
Simple 

W/0 
With 
%  dif. 

14. 

14. 

0. 

764 

762 

.013 

34.236  J 
34.232 
0.011  i 

75.065 

75.008 

0.076 

123. 

|  123. 

0. 

l 

283 
147 
110 

1  152. 

152. 

0. 

752 
516 
155 

Pinned 
Pinned 

W/0 
With 
%  dif. 

!l5 

15 

0 

.256 
.255 
,006 

25.364 

25.359 

0.020 

75.667 

75.607 

0.078 

'  149 
:.  149 

!   0 

1 

.820 
.583 
.158 

278 

277 

0 

.160 
.296 
.310 

Pinned 

i 

Simple 

W/0 
With 
%  dif. 

9 

9 
0 

.741 
.740 
.001 

24.417 

24.413 

0.016 

72.502 

72.415 

0.079 

\   122 

122 

0 

.956 
.819 
.112 

152 

152 

0 

.412 
.177 
.159 

Comparison  of  Models 

Both  the  lumped-parameter  and  finite-element-models 
lend  themselves  favorably  to  computer  analysis.   The  param- 
eters needed  for  analysis  are  easily  derived  and  fit  conven- 
iently into  a  matrix  form  which  is  easily  programmed.   Each 
model  has  one  feature  which  may  cause  some  difficulty.   These 
are    the  application  of  end  conditions  in  the  lumped- 
parameter-model   and  the  corner  transformation  in  the  finite- 
element-model.   Once  understood,  however,  both  of  these  can 
be  programmed. 

The  inertial  and  stiffness  matrices  resulting  from  these 
models  are  quite  different.   The  matrices  from  the  finite- 
element-model  are  sparse.   The  matrices  from  the  lumped- 
parameter-model  do  not  have  this  property.   The  eigenvalue 
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solution  will  therefore  be  easier  and  less  time  consuming  for 
the  finite-element-model  since  the  matrices  are  already  par- 
tially diagonalized. 

The  order  of  the  resulting  matrices  is  also  different. 
The  order  of  the  finite-element-matrices  is  2N+3  where  N  is 
the  number  of  elements.   The  order  of  the  lumped-parameter- 
matrices  is  N+l.   The  finite-element-model  using  12  elements 
or  27  degrees  of  freedom  yields  much  better  solutions  than 
the  lumped-parameter-model  with  3  8  segments  or  39  degrees  of 
freedom.   This  significant  difference  in  the  order  of  the 
matrices  will  save  computation  time  also. 
Conclusions 

1.  Both  models  are  conveniently  treated  by  digital  computer 
methods  and  could  be  used  for  the  analysis  of  frames. 

2.  Both  models  can  include  rotatory  inertia. 

3.  Both  models  provide  a  basis  for  finding  dynamic  response 
of  the  structure  to  arbitrary  time-varying  external  loads. 

4.  A  comparison  of  the  results  of  an  analysis  of  a  gable 
bent  using  the  models  shows  that  the  finite-element-model 
gives  much  better  results.   It  requires  12  fewer  degrees  of 
freedom  and  its  average  difference  from  the  standard  is  one- 
tenth  that  of  the  lumped-parameter-model. 

5.  The  transfer  matrix  method  still  gives  good  results  for 
frequencies  higher  than  the  fifth,  but  not  mode  shapes.   The 
finite-element-model  gives  good  agreement  at  these  higher 
frequencies,  and  its  mode  shapes  shooild'  fee  ;  accurate  since  the 
particular  round-off  error  that  affects  the  transfer  matrix 
method  does  not  pertain  to  the  finite-element-model. 
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APPENDIX  A 


PRINCIPAL  MODE  COMPARISON  DATA 
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TRANSFER  MATRIX  METHOD 
X  COQRDI INI 


0.0 
0.0 

180.00000 

360.00000 

540.00000 

770.00000 

720.00000 

720.00000 
CLAMPED  CLAMPEO  ENO  CONDITION 
ROTATORY  INERTIA  IS  NOT  INCLUDED 


V  COORD!  INI 

60.00000 
120.00000 
156.00000 
192.00000 
156.00000 
120.(10000 

60.00000 
0.0 


MODULUS  OF  FLASTICITYIPSI I 
10000000.00000 
^0000000.00000 
50000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 


SPECIFIC  WTILR/CU  IN) 
0.28300 
0.78300 
0.28300 
0.28300 
0.28300 
0. 28300 
0.2B300 
0.28300 


APFAISQ     |M) 

5  .88000 

S.  80000 
29.40000 
29. 40000 
29.40000 
29.40000 

4.62000 

4.62000 


MOMFNTI I w  4| 

82.50000 

»?.5OO00 

625.T0O00 

675.00000 

62"=. 00000 

625.00000 

30.30000 

10.70000 


NATURAL  FREOUENCIESIRAD/SECI 


NORMALIZED  X  DISPLACEMENTS 


0.0 
0.0 
0.0 
0.0 
0.0 


2 

3 

4 

5 

6 

7 

8 

9 

■0, 

,0365 

-0. 

1476 

-0, 

.2313 

-0. 

,2578 

-0. 

.3016 

-0. 

,3680 

-0. 

,205! 

-0, 

,0000 

0, 

.1980 

0. 

4195 

0, 

.3988 

0. 

.3440 

0. 

.3467 

0. 

,2684 

0, 

,0932 

0, 

,0000 

0, 

.1604 

0, 

1948 

0, 

.0614 

0, 

,1813 

0, 

.0422 

0, 

,1678 

0, 

,1616 

-0, 

,0000 

0. 

.0587 

-0. 

,1077 

-0. 

,21  29 

0. 

,0127 

0, 

.2284 

0. 

,1330 

-0. 

,0551 

-0, 

.0000 

0. 

.6216 

0. 

,0226 

-0, 

.0149 

0. 

,0352 

0. 

.0435 

0, 

,0479 

0. 

,7535 

-0, 

.0014 

NORMALIZED  Y  DISPLACEMENTS 


1 


0.0 
0.0 
0.0 
0.0 
0.0 


2 

0.0000 
-0.0000 
-0.0000 
-0.0000 
-0.0000 


3 

0.0000 
-0.0000 
-0.0000 

0.0000 
-0.0000 


0.4183 
0.1038 
0.6670 
0.5262 
0.1875 


5 

0.5509 
0.3777 
0.0674 
-0.6017 
-0.0633 


6 
0.3317 
0.3912 

-0.6279 
0.4770 

-0.0219 


-0.0000 
0.0000 
-0.0000 
-0.0000 
-0.0000 


-0.0000 
0.0000 
-0.0000 
-0.0000 
-0.0000 


-0.0000 
O.0000 
-0.0000 
-0.0000 
-0.0000 


LUMPFD  PARAMETER  METHOD 
X  COORD! INI 

0.0  120.00000 

360.00000  192.00000 

720.00000  120.00000 

720.00000  0.0 

CLAMPED    CLAMPEO    END    CONDITION 
ROTATORY     INERTIA     IS    NOT     INCLUDED 

NATUPAL  FRE0UENC1 ESIRAD/SECI 


NUMBER    OF     SECTIONS 

6.00000 

13.00000 

13.00000 

6.00000 


3 
79.366 


SPECIFIC    HTIL8/CII     IN) 
0.28300 
0.28100 
0.28300 
0. 28300 


ARFMSO     IN) 

5  .88000 
29.40000 
29.40000 

4.67000 


MOMFNTI IN     4 ) 

87.50000 

675.00000 

625.00000 

10.30000 


NORMALIZED    X    DISPLACEMENTS 

1 


0.0 
0.0 
0.0 
0.0 
0.0 


2 

3 

4 

5 

6 

7 

a 

q 

0, 

.0352 

0. 

,1452 

0, 

,2294 

0, 

,  2558 

0. 

.2999 

0. 

,3664 

0, 

.7042 

-0, 

,0000 

0, 

,1969 

0. 

.4192 

0, 

.4000 

0, 

,  3444 

0. 

,  3483 

0, 

,  7696 

0. 

,0934 

0, 

.0000 

0, 

,1613 

-0. 

1984 

-0, 

.0653 

-0, 

.  1845 

-0. 

,0455 

-0. 

,1706 

-0, 

,1671 

0, 

,0000 

0. 

,0600 

-0, 

1054 

-0. 

,2119 

0. 

.0134 

0. 

.2784 

0. 

,1323 

-0, 

.0560 

-0, 

,0000 

0, 

.6208 

0. 

,0247 

-0, 

.0160 

0, 

.0379 

0, 

.0464 

0. 

,0511 

0. 

,7497 

0. 

,0000 

NORMALIZED  Y  DISPLACEMENTS 


I 


0.0 
0.0 
0.0 
0.0 
0.0 


0.0 
0.0 
0.0 
0.0 
0.0 


FINITE  FLFMFNT  METHOD 

n8!8 

0.0 

360.00000 

720.00000 

720.00000 
CLAPPED    CLAMPED    ENO    C.nnolTIPN 
ROTATORY     INERTIA     IS    NOT     INCLUOEO 


7 

4 

6           6 

7           8 

Q 

0.0       -0, 

,420" 

-0.55H    -0.3376 

0.0000    -0.0000 

-0 

.0000 

0.0        0, 

,0964 

0.7741     0.1936 

-0.0000    -0.0000 

-0 

.0000 

0.0       -0, 

.6654 

-0.0693     0.6255 

0.0000     0.0000 

0 

.0000 

0.0        0, 

,5324 

-0.5943     0.4804 

-0.0003     0.0000 

0 

.0000 

0.0        0, 

.2034 

-0.0658    -0.0232 

0.0000     0.0000 

0 

,0000 

Y  COTROI IN) 
0.0 

NUMBER  OF  SECTIONS 

SBFCIFIC  WTILB/CU  INI 

ARFAISO  IN) 

2.00000 

0.28700 

5.88000 

120.00000 

4.00000 

0.28700 

29.40000 

192.00000 

4.00000 

0. 28700 

79.40000 

120.00000 

0.0 

2.00000 

0.23700 

4.62000 

MO"FNT( IN    4) 

"7.50000 

675.00000 

675.00000 

70. 70000 


N1TURAL     FREOUENCl ESIRAD/SFC I 


'ALIZEO    X    DISPLACEMENTS 


1 

0.0 

-0.0366 

-0, 

.  1476 

-0.2313 

-0, 

.2578 

-0, 

.3016 

-0. 7680 

-0.2061 

0.0 

2 

0.0 

0.1980 

0, 

,4196 

0. 7988 

0, 

.3440 

0 

.3468 

0.2685 

0.0912 

0.0 

3 

0.0 

-0.  l<-04 

-0, 

,  1946 

-0.0612 

-0 

.1812 

-0, 

.0421 

-0.1677 

-0. 1615 

0.0 

4 

0.0 

-0.0587 

0. 

,  1077 

0.2130 

-0. 

.0126 

-0, 

.  2284 

-0. 1130 

0.0551 

0.0 

6 

0.0 

-0.6213 

-0. 

,0222 

0.0183 

-0, 

,0348 

-0 

.0431 

-0.0475 

-0.7*79 

0.0 

PMAl 

I/EO  Y 

DI SPl ACEMENTS 

1 

2 

3 

4 

5 

6 

7 

Q 

1 

0.0 

0.0 

0, 

,0 

0.41 83 

0. 

6509 

0. 

,1117 

-0.0000 

0.0 

0.0 

7 

O.J 

0.0 

0. 

0 

0. 1036 

0. 

3776 

0, 

,3912 

0.0000 

0.0 

0.0 

3 

0.0 

0.0 

0, 

0 

-0.6671 

-0, 

0674 

0. 

,6779 

-0.0000 

0.0 

0.0 

4 

0.0 

0.0 

0. 

0 

-0.5267 

0. 

6016 

-0, 

,4771 

-0.0000 

0.0 

o.o 

S 

0.0 

0.0 

0. 

0 

-0. 1872 

0. 

06  32 

0, 

,0217 

-0.0000 

0.0 

0.0 
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TRANSFER    MATRIX    METHOD 

COQRQ(IN)                         YCnnROlINI        MODULUS    OF    FL  AST  IC I TY  (  PS  I  )        SPECIFIC    WTtlB/CII     TV  I        ARFUKQ     fN|  momcnti)"'    u  | 

'ooooooo.  ooooo  o.'flioo           *.8<>nio  "'.'win 

30000000.00000  0.21300               8.08000  H'.tnO"1 

'0100000.00000  0. 28*00  2",  40000  "'.^OTl 

'1000000.00000  0.20300  73.40000  675.1001" 

'0000000.00000  0.28300  ?1.«0nnn  626.00010 

'0000000.00000  0.28'00  '1.40000  fiyK.DOOOri 

30000000.00000  0.^8300               4.67000  30.30010 

30000000.00000  0.73300               4.157000  30. '0111 


0.0  60. OOOOO 

0.0  120.00000 

lao. OOOOO  15*. OOOOO 

360.00000  192.00000 

540.00000  156.00000 

720.00000  120.00000 

720.00000  60.00000 

720.00000  0.0 
CLAMPED  SIMPLE  ENO  CONDITION 

ROTATORY  INERTIA  IS  NOT  INCLUDED 

NATURAL  FREQUENCIESIRAD/SEC) 

1  2 

1  14.765  34.238 

NORMALIZED  X  DISPLACEMENTS 

1 

1  0.0 

2  0.0 


2 

3 

4 

5 

6 

7 

Q 

Q 

0. 

.0506 

-0. 

,  1*37 

-0. 

,7088 

-0, 

,2302 

-0. 

,2538 

-0.3067 

-0. 

,4073 

-0, 

.5108 

0. 

.  1953 

0. 

,4069 

0. 

.  3781 

0. 

,3217 

0, 

,  3200 

0.7366 

0, 

.0658 

-0 

.1062 

0 

.0839 

0. 

.0960 

0, 

.0186 

0. 

,0800 

-0. 

,0027 

0.0641 

0, 

,  3640 

0, 

.7571 

0, 

,0108 

0. 

.0082 

-0. 

,0011 

0. 

,0153 

0, 

,0191 

0.0225 

-0. 

.3178 

-0, 

,945'= 

0, 

.0236 

-0. 

,0561 

-0. 

.1031 

-0, 

.0016 

0. 

,1011 

0.0531 

0, 

.7574 

0, 

.0480 

3  0.0 

4  0.0 

5  0.0 

NORM4LI2E0  Y  DISPLACEMENTS 

123456789 

1  0.0                     0.0000  O.OOCO  0.2754  0.3373  0.2644  0.0000  0.0000             0.0000 

2  0.0                  -0.0000  -O.OOCO  0.1438  0.4256  0.4172  -0.0000  -0.0000  -0.1OOO 

3  0.0                  -0.0000  -0.0000  0.3S66  0.0797  -0.3337  0.0000  0.0000             0.0000 

4  0.0                  -0.0000  -0.0000  0.0463  -0.0357  -0.0171  0.0000  0.0000             0.0000 

5  0.0                  -0.0000  0.0000  0.2339  -0.2734  0.2403  0.0000  0.0000             0.0100 


LUMPED    PARAMETER    METHOD 

X    COQRQIINI                            Y    COnRD(IN)                      NUMBER    OF     SECTIONS  SPECIFIC    WTIIB/CU  IN)              4RF»(S0     TM)  »n»^:TI"i    '. ) 

6.00000                                    0.28300                                    5. ooooo  8J.5O000 

13.00000                                       0.28300                                    79.40000  675.01011 

13.00000                                       0.28300                                    2°. 40000  675.10000 

6.00000                                       0.28300                                       4.67000  ->0. 30111 


0.0  120.00000 

360.00000  192.00000 


720.00000  120.00000 

720.00000  0.0 


CLAMPED    SIMPLE    END    CONOITION 
ROTATORY     INERTIA     IS    NOT     INCLUDED 


NATURAL  FREOUENCI ESIRAD/SEC) 

1 

1 
14. 

769 

2 

34.482 

NORMALISED  X 

DISPLACEMENTS 

1 

1  0.0 

2  0.0 

3  0.0 

4  0.0 

5  0.0 

2 

0.0*00 
0.1°43 

-0.0054 
0.0110 

-0.0213 

3 
0.  1526 
0.4067 
-0.0992 
0.0075 
0.0507 

NDRMAL  1  7.50  Y 

DIS 

PLACEMENTS 

1 

1  0.0 

2  0.0 

3  0.0 

4  0.0 

5  0.0 

2 

0.0 
0.0 
0.0 
0.0 
0.0 

3 

0.0 
0.0 
0.0 
0.0 
0.0 

5 
150.914 


4  5  6  7           8           9 

0.2079  0.2293  0.2523  0.3060     0.4070     0.5108 

0.3793  0.3272  0.3216  0.7377     0.0654  -0.1080 

-0.0210  -0.0327  0.0010  -0.0662  -0.3639  -0.7526 

-0.0022  0.0152  0.0201  0.0229  -0.3174  -0.9454 

0.0930  0.0072  -0.0904  -0.0463  -0.7657  -0.8690 


-0.2767  -0.3335  -0.2663  0.0000  0.0000  0.0000 

0.1371  0.4225  0.4195  -0.0000  -0.0000  -0.0000 

-0.3906  -0.0823  0.3363  -0.0001  -0.0000  -0.0000 

0.0488  -0.0385  -0.0140  0.0000  -0.0000  -0.0000 

-0.2117  0.2423  -0.2200  -0.0000  -0.0000  -0.0000 


F I N I T  F  ELEMENT  "ETHOD 

X  CIORD(IN)  Y  rnnpn(IN)         NU-BER  OF  SECTIONS   S'FCIFIC  WTILR/CU  IN)      A9F4|<;0  r»l        »n«l=Mr(TN  4) 

120>'00000  '-OOOOO  0.28300  5.0,000  -..mil 

4.00000  1.78300  '9.40000  676.100H 

4.00001  0.78300  79.40000  ',"=.11001 
720.00000               CO                    '.OOOOO               0.78'10               4.67100  30. '1111 

CLAMFED  SIMPLE  END  CONDITION 
R0TATO3Y  INERTIA  IS  NOT  INCLUDED 


3U. OOOOO  192.00000 

720.00000  120.00000 


NATURAL  FOF-juENCIFSIRAD/SrC) 

1  2 

1  14.764  34.236 


►  inKMALI7F0    X    DI  SPl  ACf  MFNTS 

1 

1  0.0 

2  0.0 

3  0.0 

4  0.0 

5  0.0 

'.DRMALIZEO  Y  riSPLACfMFNTS 


' 

3 

4 

6 

6 

7 

o 

Q 

-0, 

.0*0* 

-0 

.1637 

-0, 

,7038 

-0, 

.2102 

-0. 

.2638 

-0, 

.  3067 

-0 

,40'3 

-0 

.*'  10 

0, 

,  1953 

0 

.4068 

0, 

,  1731 

0. 

,3717 

0, 

,  3200 

0. 

,7366 

0 

.0667 

-0, 

.1  06? 

-0 

.0839 

-0 

.  0969 

-0, 

,01  86 

-0. 

,0800 

0. 

,0027 

-0, 

,0640 

-0, 

,  '640 

-0, 

.7671 

-0. 

,0108 

-0 

.00P2 

0. 

.0011 

-0. 

,0153 

-0. 

0  110 

-0. 

.0726 

0 

.3170 

0, 

,  9465 

0, 

.0236 

-0, 

.0664 

-0, 

,1033 

-0, 

,0016 

0, 

1013 

0, 

.0637 

0, 

.2*7? 

0, 

,84'A 

1 

2  0.0  0.0 


0.0  0.0  0.0 

0.0  0.0  0.0 


4 

6 

6 

7 

0. 

.7784 

0, 

.  1873 

0, 

,2644 

-0.3000 

0 

.0 

0, 

.1438 

0. 

,4266 

0, 

,4177 

J. 0000 

0, 

.0 

0. 

,  386* 

-0, 

,  0'97 

0. 

,3337 

-0.0000 

0, 

,  0 

0. 

0463 

0, 

,0367 

0. 

,017! 

-0.0000 

0, 

.0 

0, 

,2  345 

-0, 

,7739 

0. 

,2407 

0.0000 

0, 

.0 

0.1 
0.0  0.0  0.0  -0.386*  -C0797  513337  -COOOO  0 !  0  o! 


.0 
0.0 
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TRANSFER  MATRIX  METHOD 

x  cnpRni  INI 

0.0  60.00000 

0.0  120.00000 

180.00000  156.00000 

360.00000  192. OOOOO 

540.00000  156.00000 

720.00000  120.00000 

720.00000  60.00000 

720.00000  0.0 
PINNED  PINNED  ENO  CONDITION 

ROTATORY  INERTIA  IS  NOT  INCLUDED 


MODULUS  OF  ELASTICITYIPSI I 
30000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 
■"0000000. OOOOO 
30000000.00000 


SPECIFIC  WTUR/CU  IN)   ARFAISO  I*') 

MOMFMTI |»  1| 

0.28300 

5.88000 

«2.5oonn 

0.28300 

5.88000 

R',50000 

0.'8300 

29.40000 

625. OOOOO 

0.28300 

29.40000 

625.00003 

0.28300 

29.40000 

625. OOOOO 

0.28300 

29.40000 

625.00000 

0.28300 

4.62000 

30. 30000 

0.28300 

4.62000 

30.30000 

NATURAL  FREOUENCIESIRAD/SECI 


NORMALIZED  X  DISPLACEMENTS 


1 


0.0 
0.0 
0.0 
0.0 
0.0 


2 

3 

4 

5 

6 

7 

8 

9 

0, 

.2170 

0. 

,3622 

0. 

,3957 

0. 

,3960 

0, 

.4152 

0, 

.4290 

0, 

,2982 

-0, 

,0000 

0. 

.2132 

0. 

.2569 

0. 

,1777 

0. 

,  1257 

0, 

,0927 

-0, 

,0054 

-0. 

.0689 

0 

.0000 

0. 

,1993 

0. 

,1320 

-0. 

,0020 

0. 

,1233 

-0. 

,0153 

0, 

.1145 

0. 

,2040 

-0. 

.0000 

0, 

.0961 

-0, 

,1133 

-0, 

,2159 

0. 

,0078 

0. 

.2238 

0. 

,1289 

-0, 

.0993 

-0, 

,0000 

0. 

.6120 

-0. 

,0077 

-0, 

.0229 

-0, 

,0013 

0. 

,0061 

0, 

.0052 

0. 

.7861 

0, 

,0000 

NORMALIZED  Y  DISPLACEMENTS 


0.0 
0.0 
0.0 
0.0 
0.0 


2 

-0.0000 
-0.0000 
-0.0000 
-0.0000 
-0.0000 


-o.ooco 

-0.0000 

-0.0000 

0.0000 

0.0000 


-0. 1678 
0.3961 
0.6702 
0.5130 
0.0758 


-0.1693 
0.6559 
0.04  38 
-0.6055 
-0.0323 


6 

-0.0733 
0.4906 

-0.6493 
0.4743 
0.0046 


-0.0000 

0.0000 

-0.0000 

-0.0000 

0.0000 


-0.0000 
0.0000 
-0.0000 
-0.0000 
0.0000 


-0.0000 

0.0000 

-0.0000 

-0.0000 

0.0000 


LUMPED  PARAMETER  METHOD 
X  COORDIIN) 


Y  COORD! IN) 


0.0  120.00000 

360.00000  192.00000 

720.00000  120.00000 

720.00000  0.0 

PINNED  PINNEO  END  CONDITION 
ROTATORY  INERTIA  IS  NOT  INCLUDED 


NUMBER  OF  SECTIONS 

6.00000 

13.00000 

13.00000 

6.00000 


SDFCIFIC  WTIL8/CU 
0.28300 
0.28300 
0.2»*00 
0.28300 


ARCAI SO  IN) 

5.88000 
29.40000 
2O.40000 

4.62000 


Mr,»rMT(  IN  41 

R?.cOono 

625.00000 

625.00000 

'0.30000 


NATURAL  FREOUENCIESIRAD/SECI 


1 

1 
15.285 

2 

25.365 

3 

75 

.162 

4 
147.574 

5 
269.188 

NORMAL  1  ZED  X 

DISPLACEMENTS 

1 
2 
3 
4 
5 

1 
0.0 
0.0 
0.0 
0.0 
0.0 

2 

-0.2169 
-0.2140 
-0.2003 
-0.0991 
-0.6146 

-0.3618 

-0.2574 

-0.1327 

0.1116 

0.0081 

-0, 

-0. 

0. 

0, 

0, 

4 
.3956 
,1783 
,0013 
.2152 
.0243 

5 
-0. 3958 
-0. 1264 
-0.1238 
-0.0081 
0.0015 

6 
-0.4151 
-0.0936 
0.0148 
-0.2234 
-0.0062 

7           8 
-0.4297    -0.5937 
0.0046     0.0689 
-0.1150    -0.2048 
-0.1279     0.1023 
-0.005?    -0.7835 

0, 

0, 
0, 
0, 
0, 

Q 
,0000 
,0000 
.0000 
.0000 
.0000 

NORMALIZED  Y 

DI SPLACFMENTS 

1 
2 

4 

5 

1           2 
0.0        0.0 
0.0        0.0 
0.0        0.0 
0.0        0.0 
0.0        0.0 

cLFMfNT  METHno 

3 
0.0 
0.0 
0.0 
0.0 
0.0 

0. 
-0, 
-0. 
-0. 
-0. 

4 

.1689 
.3953 
,6698 
,51  73 
.0807 

5 

0.  1698 

-0.6550 

-0.0444 

0.5938 

0.0332 

6 

0.0733 
-0.4913 

0.6488 
-0.4775 
-0.0054 

7           8 
-0.0000    -0.0000 

0.0000     0.0000 
-0.0000    -0.0000 

0.0000     0.0000 
-0.0000    -0.0000 

-0, 
0, 

-0. 
0, 

-0. 

9 

.0000 
.0000 
,0000 
.0000 
.0000 

FINITE 

X  COORD! IN) 

Y  COORDI 

IN) 

NUMRER  OF 

SECTIONS 

SPECIFIC  WTILR/CU  IN) 

ARFA<«0  !►') 

MOMpKjTl  [►'  4) 

0 
360, 
720, 
720, 

,0 

.00000 
.00000 
.00000 

120.00000 
192.00000 
120.00000 
0.0 

2.00000 
4.00000 
4.00000 
2.00000 

0.28300 
O.'»300 
0.28300 
0.23300 

5.88000 

?i.4nooo 

2°.4O000 

4.62000 

8', 

625, 
6?5, 

30, 

.80000 
,30000 
,30000 

,30000 

PINNED 

PINNEO  END  CONDITION 

ROTATORY  INERTIA  IS  NOT  INCLUDED 

NATURAL  FREQUENCIES  I  SAD/ SEC ) 


NORMALIZED  X  DISPLACEMENTS 


0.0 
0.0 
0.0 
0.0 
O.U 


2 

3 

4 

5 

6 

7 

8 

J, 

.2170 

0. 

.  3622 

0, 

.  3957 

0, 

3960 

0, 

.4152 

0, 

,4299 

0 

,298? 

0, 

,0 

0, 

,2132 

0, 

,  2569 

0, 

.1777 

0. 

,1257 

0. 

,09?7 

-0, 

,0054 

-0, 

,0689 

0, 

,0 

0, 

.  1993 

0. 

,  1319 

-0, 

,0021 

0, 

,1232 

-0, 

,0154 

0, 

,  1  144 

0, 

,2040 

0, 

,0 

0, 

,0961 

0. 

,1133 

0, 

,2159 

-0. 

,0078 

-0. 

,2238 

-0, 

,1289 

0, 

,0903 

0, 

,0 

0, 

.6119 

0, 

,0078 

0. 

,0229 

0. 

0013 

-0, 

,0060 

-0, 

0051 

-0, 

,786? 

0, 

,n 

■<MALt/ED    Y     DI  SPL  ACCENTS 


1 


0.0 
0.0 
0.0 
0.0 
CO 


0.0 
0.0 
0.0 
0.0 
0.0 


0.0 
0.0 
0.0 


0, 

4 
.1678 

-0, 

5 
,  1693 

-0, 

6 

,0733 

0, 

7 
,0000 

0 

,0 

0. 

.  3961 

0, 

,6559 

0, 

,4906 

-0, 

.0000 

0 

.0 

n 

.6702 

0, 

,04  38 

-0, 

,6493 

0, 

.0000 

0 

,n 

0, 

.51  30 

0. 

,6054 

-0, 

,4743 

-0, 

.0000 

0, 

.0 

0, 

.0756 

0, 

,0322 

-0. 

,0046 

-0, 

,0000 

0, 

.0 

0.0 
0.0 
0.0 
0.0 
0.0 
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TRANSFER  MATRIX  METHOD 
X  COORD! IN) 


0.0 
0.0 

180.00000 

360.00000 

540.00000 

720.00000 

720.00000 

720.00000 
PINNED    SIMPLE    ENO    CONDITION 
ROTATORY     INERTIA     IS    NOT     INCLUDED 


V    nilRfUIN) 

60.00000 
120.00000 
156.00000 
192.00000 
156.00000 
120.00000 

60.00000 
0.0 


MODULUS  OF  ELAST1CITYIPSI > 
30000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 


SPECIFIC  WT(LH/r.U  INI 
0.23301 
0.28300 
0.28300 
0.28300 
0.28300 
0.28300 
0.28300 
0.28300 


ARFMSO     TN| 

5.88000 

^.88000 
29.40000 
29.40000 
29.40000 
29.40000 

4.62000 

4.6'onn 


MOMENT!  IN  4  1 

87.50010 

8->.<:oP00 

625. OOOOO 

6'*. OOOOO 

625. OOOOO 

ft'S. 31001 

30.^0000 

3l.3iooo 


NATURAL  FREQUENCIESIRAD/SECt 


3 
72.489 


NORMALIZED  X  DISPLACEMENTS 


1 


0.0 
0.0 
0.0 
0.0 
0.0 


2 

3 

4 

5 

6 

7 

8 

9 

-0 

.  1556 

-0, 

.2682 

-0. 

.3028 

-0. 

3120 

-0. 

.3270 

-0, 

.3558 

-0. 

.4092 

-0, 

.4638 

0 

.2354 

0. 

,3043 

0. 

,2398 

0. 

,  1903 

0, 

,  1668 

0, 

.0763 

-0, 

.1016 

-0. 

,2830 

0 

.1120 

0. 

,0700 

-0. 

,0105 

0, 

,0584 

-0. 

.0269 

0, 

.0469 

0, 

.3520 

0, 

.7421 

0 

.0152 

0. 

.0061 

-0. 

,0027 

0. 

0137 

0. 

.0183 

0. 

.0213 

-0, 

.3189 

-0, 

.9452 

0, 

.0377 

-0. 

.0589 

-0, 

.1035 

-0. 

0051 

0. 

,0961 

0, 

.0487 

0, 

.2577 

0, 

.8532 

NORMALIZED  Y  DISPLACEMENTS 


LUMPFO  PARAMETER  METHOD 
X  COORD! IN) 

0.0 

360.00000 

720.00000 

720.00000  0 

PINNED  SIMPLE  END  CONDITION 
ROTATORY  INERTIA  IS  NOT  INCLUDED 


2 

0.0000 
0.0000 
0.0000 
0.0000 
0.0000 

3 

0.0000     0, 
-0.0000     0. 
-0.0000     0. 
-0.0000     0. 

0.0000     0. 

4 
.1732 
,3227 
.4026 
,0444 
,2230 

5 

0.2190     0. 

0.5701     0, 

0.0578    -0, 

-0.0379    -0, 

-0.2689     0. 

6 

,1439 
4527 
,3690 
,0151 
,2375 

ooooo 

CCOOO 

cr  cocoo 

ooooo 

ooooc 

ocooo 

ooooo 

*-coooo 

ooooo 

ooooo 
1 

0 
-0 
0, 
0 
0, 

9 
.0000 
,0000 
.0000 
.0000 
.0000 

IETHOD 

Y  COORD!  IN) 

120.00000 
192.00000 
120.00000 

NUMBER  OF  SECTIONS 

SPECIFIC  WT(l B/CU  IN) 

AOFAISQ  IN) 

1 
1 

6.00000 
I  3.00000 
13.00000 

6.00000 

0.28300 
0.28300 
0.28300 
0.28300 

5.08000 
29.40000 
29.40000 

4.62000 

*n"FNTI  IN  4) 

82.50000 

675.00000 

625.00000 

30.30000 


NATURAL  FREQUENC IESIRAD/SEC) 


2 

24.419 


123.511 


NORMALIZED  X  DISPLACEMENTS 


0.0 
0.0 
0.0 
0.0 
0.0 


0. 

2 
,1555 

0. 

3 
.2679 

0, 

4 
3026 

0. 

5 
,3118 

0. 

6 

,  3?68 

0. 

7 
.3557 

0. 

,2361 

-0. 

,3049 

-0. 

2406 

-0. 

,  19U 

-0, 

1680 

-0. 

,0773 

0, 

,  1140 

0, 

,0714 

-0. 

,0101 

0, 

,0597 

-0. 

,0268 

0, 

.0479 

0. 

.0156 

0. 

,0054 

-0. 

,0040 

0. 

,0135 

0. 

,0193 

0, 

.0217 

0, 

.0344 

0. 

.0531 

0. 

,0933 

0. 

,0056 

-0. 

,0355 

-0, 

.0419 

0.4092 

0. 1007 

0.3512 

-0.3186 

-0.265° 


0.463" 

0.2822 

0.7370 

-0.9451 

-0.873'' 


NORMALIZED  Y  DISPLACEMENTS 


1 


0.0 
0.0 
0.0 
0.0 
0.0 


FINITE  ELEMENT  METHOD 

0.0 

360.00000 

720.00000 

720.00000 
PINNED    SIMPLE    ENO    CONDITION 
RCTATCRY     INERTIA     IS    NOT     INCLUDED 


3 

4 

5 

6 

7           8 

9 

0.0       -0. 

,1740 

-0.2195 

-0. 1444 

0.0000     0.0000 

0, 

.0000 

0.0       -0. 

3216 

-0.5691 

-0.4534 

0.0000     0.0000 

0, 

.0000 

0.0        0. 

,4073 

0.0587 

-0.3734 

0.0000     0.0000 

0, 

.0000 

0.0        0, 

,0470 

-0.0409 

-0.0118 

0.0000    -0.0000 

-0 

.0000 

0.0       -0. 

,2006 

0. 2376 

-0.2178 

0.0000     0.0000 

-0, 

.0000 

y  coren( in) 

Nl 

l»8FR  OF 

SECTI ONS 

SPFC1FIC  WTUR/PI  IN) 

49 F  SI  SO  IN) 

2.00000 

0.2R300 

s .08000 

120.00000 

4.0000 

0.28300 

2°.4000O 

192.00000 

4.0000 

0.28300 

29.40000 

120.00000 

2.0000 

0. 28300 

4.67000 

0.0 

FOMENT!   |M     4) 
07.50100 

6?":. 00000 

6'5. 30010 
30.30001 


NATURAL  FREOUENC I ESIRAQ/SEC) 


2 

24. <-17 


5 
152.412 


N.lfMAIIZEO    X    DISPLACEMENTS 


2 

3 

4 

5 

6 

7 

8 

9 

0. 

,  1556 

0, 

.2682 

0, 

,3028 

0. 

,  3120 

0. 

.3270 

0, 

.  3558 

0, 

.409? 

0, 

.4630, 

0, 

.2354 

0, 

,  3043 

0. 

2398 

0. 

1903 

0. 

,  1668 

0, 

.0763 

-0, 

,1017 

-0 

.2830 

-0. 

.1119 

-0, 

.0699 

0. 

.0106 

-0. 

,0584 

0. 

.0270 

-0, 

,0468 

-0. 

,38->0 

-0, 

.7421 

-0, 

.0152 

-0, 

,0061 

0. 

,0027 

-0, 

,0137 

-0, 

.0182 

-0 

.071  3 

0, 

.3189 

0 

.9457 

-0, 

.0374 

0. 

,059Q 

0. 

,1037 

0, 

0051 

-0, 

,0964 

-0, 

.0488 

-0. 

,2574 

-0, 

.8528 

NORMALIZED    Y    DISPLACEMENTS 


0.0 
0.0 
0.0 
0.0 
0.0 


0.0 
0.0 
0.0 
0.0 
0.0 


4 

5 

6 

7 

0, 

,1732 

-0. 

,2190 

-0, 

.  1439 

0, 

,0000 

0, 

,0 

0. 

3227 

0, 

,6702 

0. 

,4527 

0, 

,0000 

0, 

.0 

0, 

,4025 

-0, 

,0578 

0. 

,3690 

-0, 

,0000 

0. 

.0 

0. 

,0444 

0, 

,0378 

0. 

.0152 

-0. 

,0000 

0, 

,0 

0. 

2236 

0. 

,7694 

-0. 

,2379 

-0. 

,0000 

0, 

,0 

0.0 
0.0 
0.0 
0.0 
0.0 
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FINITE  ELEMENT  METHOD 

X    COORQIINI                           Y    COORDdN)                      NUMBER    OF    SECTIOMS  SPECIFIC    HTILB/CH  INI             ARFAISO    INI  MOMFNTMN    <.  I 

2.00000              0.28300              5.8BOO0  R?.5o™« 

4.  00000              0.26300             25.40000  6?5.1010n 

4.00000              0.28300             29.40000  625 . 00100 

2.00000              0.28300              4.62000  30. 'nm 


0.0  120.00000 

360.00000  192.00000 


720.00000  120.00000 

720.00000  0.0 


CLAMPED  CLAMPEO  ENO  CONDITION 
ROTATORY  INERTIA  IS  INCLUDED 

NATURAL  FREQUENCIESIRAO/SECI 

1  2 

1  25.748  34 

NORMALIZE!?  X  DISPLACEMENTS 

1 

1  0.0 

2  0.0 

3  0.0 

4  0.0 

5  0.0 

NORMALIZED  Y  DISPLACEMENTS 

12345678C 

1  0.0  0.0  0.0  0.4183  0.5510  0.3318  -0.0000  0.0  0.0 

2  0.0  0.0  0.0  0.1034  0.3774  0.3111  0.0000  0.0  0.0 

3  0.0  0.0  0.0  0.6670  0.0675  -0.6278  0.0000  0.0  0.0 

4  0.0  0.0  0.0  0.5264  -0.6015  0.4771  0.0000  0.0  0.0 

5  0.0  0.0  0.0  0.1878  -0.0633  -0.0219  0.0000  0.0  0.0 


2 

3 

4 

5 

6 

7 

8 

0 

,0364 

-0. 

,  1475 

-0, 

,2312 

-0. 

.2577 

-0. 

,3016 

-0, 

,3679 

-0, 

,2051 

0, 

.0 

0. 

,  1980 

0. 

,4196 

0, 

.3989 

0. 

.3441 

0. 

,  3469 

0. 

2686 

0, 

.093? 

0, 

,0 

0, 

.1605 

0. 

,  1948 

0. 

,0614 

0, 

.1813 

0, 

,0423 

0, 

,1676 

0, 

,161*. 

0, 

,0 

0, 

,0587 

-0. 

,1076 

-0, 

,2129 

0. 

,0127 

0. 

,2284 

0, 

,1330 

-0 

.0552 

0, 

,  o 

0. 

.6214 

0, 

,0225 

-p, 

,0151 

0, 

.0351 

0, 

,0434 

0, 

,0478 

0, 

,753ft 

<"> , 

.6 

FINITE  ELEMENT  METHOD 

X    COOROIINI  Y    COORQ(IN)  NUMBER    OF    SECTIONS       SPECIFIC    WTILB/CIJ    INI  AREAtSO    INI  MtJMFNTIIN    4| 

2.00000  0.28300  5.08000  82.S010O 


4.00000  0.28300  29.40000  ft?5.non"P 

4.00000  0.28300  29.40000  625.00010 

2.00000  0.28300  4.62000  lO.'OOin 


0.0  120.00000 

360.00000  192.00000 

720.00000  120.00000 

720.00000  0.0 
CLAMFEC  SIMPLE  ENC  CONDITION 
ROTATORY  INERTIA  IS  INCLUDED 

NATURAL  FPEQUFNCI ESIRAC/SEC I 

12  3  4  5 

1  14.762  34.232  75.008  123.147  152.516 

NORMALIZEH  X  DISPLACEMENTS 

1 

1  0.0 

2  O.C 

3  0.0 

4  0.0 

5  CO 

NCRMALIZEC  Y  DISPLACEMENTS 

12345678C 

1  0.0  0.0  0.0  0.2754  0.3823  0.2644  -0.0000  0.0  0.0 

2  0.0  0.0  0.0  0.1437  0.4255  0.4171  0.0000  0.0  0.0 

3  CO  0.0  0.0  -0.3866  -0.0798  0.3336  -0.0000  0.0  0.0 

4  0.0  0.0  CO  -0.0463  0.0358  0.0171  -0.0000  0.0  0.0 

5  0.0  0.0  0.0  -0.2340  0.2733  -0.2403  -0.0000  0.0  0.0 


2 

3 

4 

5 

6 

7 

8 

9 

-0, 

,0506 

-0, 

,  1537 

-0, 

,2088 

-0, 

,2302 

-0. 

,2538 

-0, 

,  3066 

-0, 

,4073 

-0. 

,5108 

0, 

.1953 

0. 

,4069 

0. 

,3781 

0, 

,3218 

0, 

3201 

0, 

,2367 

0, 

,0658 

-0, 

.1062 

-0, 

.0839 

-0. 

,0960 

-C. 

,0187 

-0. 

,0800 

0, 

,0027 

-0, 

.0641 

-0. 

.3640 

-0, 

.7572 

-0, 

.0108 

-0, 

.0082 

0. 

,0011 

-0, 

,0153 

-0. 

,0191 

-0, 

.0225 

0, 

.3177 

0, 

,9456 

-0, 

.0236 

0, 

,0563 

0. 

,1031 

0, 

,0016 

-0. 

.1011 

-0, 

.0530 

-0, 

.2573 

-0, 

,8481 
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FINITE  ELEMENT  METHOD 


x  cagnpi ini 

y 

COORD! INI 

DUMBER  OF  SECTIONS 

SPECIFIC  WT(L8/CU  INI 

AREAI SO  INI 

MOMFNTI  TN  4 

2.00000 

0.28300 

5.88000 

fl'.^nono 

0.0 

120.00000 

4.00000 

0.28300 

29.40000 

6?5. ooooo 

360.00000 

192.00000 

4.00000 

0. 28^00 

29.40000 

62c. 30000 

720.00000 

120.00000 

2.00000 

0.28300 

4.62000 

3n.iop<"n 

720.00000 

0.0 

PINNED  PINNED  END  CONDITION 

ROTATORY  INERTIA  IS 

INCLUDED 

NATURAL  FRE0UENC1ESIRAD/SECI 

1 
1           15.255 

2 
25.15<5 

1 
75 

.607 

* 

149, 

,583 

5 
277. 296 

NORMALIZED  X  DISPLACEMENTS 

123456789 

1  0.0                     0.2170  0.3621  0.3957  0.3960  0.4152  0.4299  0.2982  0.0 

2  0.0                   0.2133  0.2570  0.1777  0.1258  0.0927  -0.0054  -0.0688  0.0 

3  0.0       -0.1994  -0.1320  0.0020  -0.1233  0.0154  -0.1145  -0.2040  0.0 

4  0.0       -0.0962  0.1132  0.2158  -0.0078  -0.2238  -0.1289  0.0993  0.0 

5  0.0       -0.6123  O.C077  0.0229  0.0013  -0.0061  -0.0052  -0.7859  0.0 

NORMALIZED  Y  DISPLACEMENTS 


1 

2 

3 

4 

5 

6 

7 

8 

1 

0.0 

0.0 

0.0 

-0.1678 

-0.1693 

-0.0733 

0.0000 

0.0 

0.0 

2 

0.0 

0.0 

0.0 

0.3961 

0.6559 

0.4906 

-o.oooo 

0.0 

0.0 

3 

0.0 

0.0 

0.0 

-0.6702 

-0.04  38 

0.6493 

-0.0000 

0.0 

0.0 

4 

0.0 

0.0 

0.0 

-0.5131 

0.6053 

-0.4743 

-0.0000 

0.0 

0.0 

5 

0.0 

0.0 

0.0 

-0.0758 

0.0322 

-0.0045 

-0.0000 

0.0 

0.0 

FINITE  ELEMENT  METHOO 

X  C00R0I1NI  Y  COCRDIINI         NUMBER  OF  SECTIONS   SPFCIFIC  WT(LB/CU  INI     ARFAISQ  INI        Mr,MfNT(f  41 


0.0  120.00JOO 

360.00000  192.00000 

720.00000  12c. OOOOO 

720.00000  0.0 

PINNEO  SIMPLE  END  CONDITION 

ROTATORY  INERTIA  IS  INCLUDED 

NATUCAL  FREOUtNCI ESIRAD/SEC) 


2.00000  0.28300  5.88000  8?. 50000 

4.00000  0.28300  29.40000  6->5.0i">010 

4.00000  0.28300  29.40000  625.00000 

2.00000  0.28100  4.62O00  30.10000 


3  4  5 

72.445  122.819  152.177 


2 

i 

4 

5 

6 

7 

8 

q 

0, 

,  155t> 

0, 

,2682 

0. 

.3028 

0. 

3120 

0, 

,  3270 

0, 

,3553 

0, 

,4092 

0. 

,4638 

0, 

.2354 

0, 

,3044 

0, 

.2399 

0, 

1904 

0. 

1669 

0, 

,0763 

-0, 

,1016 

-(J, 

.2830 

-0. 

,1120 

-0. 

,C700 

0, 

,0105 

-0, 

,0584 

0. 

0269 

-0 

,0469 

-0, 

,3520 

-0. 

.7421 

-0, 

,C152 

-0, 

0C61 

0, 

,0028 

-0, 

0137 

-0. 

,0183 

-0, 

.0213 

0. 

,3188 

0, 

,94>;? 

-0, 

.0378 

0. 

,C565 

0, 

,1035 

0, 

00  5  1 

-0. 

0961 

-0, 

.0436 

-0. 

,2575 

-0 

.8533 

normalized  x  displacements 

i 

1  o.o 

2  o.o 

3  0.0 
0.0 

5      0.0 

NCRMALIZEC  Y  DISPLACEMENTS 

12345678  = 

1  0.0  0.0  0.0  -0.1732  -0.2190  -0.1439  0.0000  0.0  0.0 

2  0.0  0.0  0.0  0.3227  0.5701  0.4527  0.0000  0.0  0.0 

3  O.C  0.0  0.0  -0.4026  -0.0578  0.36°0  -0.0000  0.0  0.0 

4  0.0  0.0  0.0  -0.0445  0.0379  0.0151  -0.3000  0.0  0.0 
*              0.0  0.0  0.0  -0.2230  0.2688  -0.2375  -0.0000  0.0  0.0 
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APPENDIX  B 

LUMPED-PARAMETER  PROGRAM 
General 

Program  LUMPED  is  a  double  precision  IBM  FORTRAN  IV 

7 
language   digital  computer  program  designed  to  find  the 

principal  vibrational  modes  of  straight  sided  plane  frames 
of  uniform  cross  section.   The  program  is  written  for  the 
IBM  360/67  facility  installed  at  the  Naval  Postgraduate 
School.   The  in-plane  undamped  natural  frequencies,  as  well 
as  normalized  displacements,  may  be  found  through  its  use. 
The  displacements  are  calculated  for  the  mid-point  and  ends 
of  each  straight  member.   The  analysis  is  performed  by  the 
characterization  of  the  structure  as  a  chain  of  pin-connec- 
ted rigid  bars  with  both  mass  and  flexibility  lumped  at  the 
hinges.   The  flexibility  is  represented  by  torsional  springs. 
The  model  neglects  axial  and  shear  deformation,  but  rotatory 
inertia  can  be  included  if  the  user  desires. 

The  program  is  dimensioned  so  that  a  structure  with 
up  to  nine  straight  members  may  be  analyzed.   These  members 
may  be  divided  into  a  maximum  of  79  segments  (total) .   The 
program  with  this  dimensioning  requires  total  storage  of 
400K  bytes. 

The  more  severely  constrained  end  of  the  structure  is 
designated  the  left  end.   This  end  must  be  either  fixed  or 
pinned.   The  right  hand  end  is  allowed  any  of  the  four 
common  boundary  conditions:   fixed,  pinned,  simply  supported, 
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or  free.   The  structure  may  not  have  any  intermediate  sup- 
ports or  externally  applied  forces. 
Program  Structure 

The  program  consists  of  a  main  body  and  two  subroutines. 
If  double  precision  sine,  cosine,  arctangent,  absolute  value, 
and  square  root  routines  are  not  available  in  the  FORTRAN 
library,  they  must  be  provided  by  other  means  since  they  are 
required  by  the  program.   The  main  program  contains  comment 
cards  explaining  each  major  section's  function,  and  a  large 
leading  section  giving  input  format  information. 

When  the  end  conditions  are  applied,  the  order  of  the 
matrices  is  reduced.   Rather  than  rearrange  the  storage  of 
these  matrices,  the  rows  and  columns  corresponding  to  the 
no-longer-needed  equations  and  rotations  are  replaced  with 
zeroes.   In  the  mass  matrix,  each  diagonal  element  corre- 
sponding to  a  deleted  row  and  column  is  replaced  with  unity. 
The  eigenvalue  problem  may  now  be  undertaken.   The  final 
spectral  matrix  will  contain  zeroes  for  the  eigenvalues  of 
the  deleted  rows. 

The  subroutine  JACVAT  obtains  a  solution  for  the  real 
eigenvalue  problem  of  a  symmetric  matrix.   The  subroutine 
is  based  on  Jacobi's  diagonalization  by  successive  rota- 

g 

tions.    The  subroutine  will  take  a  matrix  up  to  the  maximum 
size  of  160  by  160.   Care  must  be  taken  to  insure  that  the 
entering  matrix  is  symmetric.   If  any  of  the  off-diagonal 
elements  are  non-symmetric,  the  subroutine  prints  the  loca- 
tion of  these  elements,  and  then  makes  them  symmetric  by 
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averaging  them.   This  output  will  be  printed  in  the  midst 
of  the  regular  output,  and  is  therefore  not  desirable. 

The  subroutine  MAPRIN  is  the  output  subroutine.   This 
subroutine  will  print  a  matrix,  column  by  column.   It  num- 
bers the  rows  and  columns  and  allows  a  single  line  heading 
of  32  spaces.   The  heading  must  be  in  a  DATA  statement  in 
the  calling  program.   The  subroutine  has  six  different 
formats  which  may  be  used. 
Input 

For  program  input,  the  following  things  are  needed: 

1.  Coordinates  of  the  ends  and  the  corners  of  the  struc- 
ture.  A  Cartesian  coordinate  system  is  used,  with  the 
positive  x  direction  to  the  right  and  the  positive  y  direc- 
tion up.   All  coordinates  are  in  inches,  and  start  with 
the  more  severely  constrained  (left)  end. 

2.  Number  of  segments  into  which  each  member  is  to  be 
divided. 

3.  Properties  of  each  member. 

a.  specific  weight,  lb/cu.in. 

b.  area,  sq.in. 

c.  modulus  of  elasticity,  psi. 

d.  second  moment  of  the  cross  sectional  area  with 
respect  to  the  z  axis,  in. 4 

4.  End  conditions. 

5.  Angle  at  which  right  hand  end  condition  is  applied, 
from  +90°  to  -90°.   This  angle  is  of  concern  only  when  the 
right  end  is  simply  supported.   For  all  other  cases,  it  is 
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taken  as  zero.   The  angle  is  the  acute  angle  between  the  y 
axis  and  the  line  of  action  of  the  restraining  force.   It 
is  measured  from  the  y  axis  and  is  positive  counterclock- 
wise . 

The  FORMAT  for  input  data,  as  well  as  the  make-up  of 
the  data  deck,  are  listed  in  the  comment  section  at  the 
start  of  the  main  program. 
Output 

The  program  will  automatically  print  for  each  problem 
and  set  of  end  conditions  within  a  problem  the  following: 

1.  Coordinates  of  the  structure  and  number  of  segments 
between  each  set  of  coordinates. 

2.  The  properties  of  each  member.   These  properties  are 
the  same  as  those  in  the  input  section. 

3.  End  conditions. 

4.  Whether  or  not  rotatory  inertia  has  been  included. 

The  user  has  control  over  the  rest  of  the  output. 
The  user  specifies  the  number  of  natural  frequencies  de- 
sired with  card  number  220  in  the  program  listing.   The 
total  number  of  frequencies  requested  may  not  exceed  the 
number  of  degrees  of  freedom  of  the  model  minus  the  number 
of  constraints.   The  specified  number  of  frequencies  will 
be  printed  in  ascending  order.   The  mode  shapes  correspond- 
ing to  the  frequencies  are  printed  in  two  arrays.   The  u 
components  are  printed  first  and  the  v  components  second. 
Each  row  of  these  arrays  corresponds  to  a  particular 
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frequency.  The  frequencies  increase  going  down  the  columns 
of  the  arrays.  Each  mode  shape  is  normalized  such  that  the 
square  root  of  the  sum  of  the  squares  is  unity. 
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LUMPED    PARAMETER    METHOD 


0.0 

360.00000 

720.00000 

720.00000 
CLAMPEO    CLAMPED    ENO    CONDITION 
ROTATORY     INERTIA     IS    NOT     INCLUDED 

NATURAL     FREQUENCIESIRAD/SECI 


Y  COPRDI IN) 

120.00000 
192.00000 
120.00000 
0.0 


NUMRFR  OF  SECTIONS 

6.00000 
13.00000 
13.00000 

6.00000 


SPECIFIC  WTILB/CU  INI 
0.28300 
0.28300 
0.28300 
0.28300 


4RFAIS0  IN) 

4.62000 
'9.40000 
29.40000 

5.88000 


MTMCMTI 1M  4 ) 

30. '0000 

625.00000 

625.30000 

82.50000 


1 

1 

25.921 

2 

34.929 

3 
79. 

366 

4 
148.262 

5 
271, 

.688 

NORMAL) 

iZED  X 

DISPLACEMENTS 

1 
2 

3 
4 
5 

1 

0.0 
0.0 
0.0 
0.0 
0.0 

2 

0.2042 
0.0934 
-0.1623 
0.0560 
0.749  2 

3 

0.3664 

0.2696 

-0. 1706 

-0. 1323 

0.0511 

0. 

0. 
-0. 
-0. 

0. 

4 
,2999 
.3483 
.0455 
,2284 
,0464 

0. 

0. 
-0. 
-0, 

0, 

5 

,2558 
,3444 
,1845 
,0134 
,0379 

6 
0.2294 
0.4000 

-0.0653 
0.2119 

-0.0160 

0 

0, 

-0. 

0, 

0, 

7 
,1452 
.4192 
,1984 
.1054 
,0247 

0, 

0, 

-0, 

-0, 

0. 

8 
.0352 
.1969 

.1613 
,0600 
,6208 

0, 
0 

-0, 
0, 

-0, 

9 
.0000 
,0000 
,0000 
,0000 
.0000 

NORMALIZED  Y 

DISPLACEMENTS 

1 
2 

3 
4 
5 

1 

0.0 
0.0 
0.0 
0.0 
0.0 

2 
0.0 
0.0 
0.0 
0.0 
0.0 

3 
0.0 
0.0 
0.0 
0.0 
0.0 

0. 
-0. 
-0, 

0. 

0, 

4 
,3326 
,3936 
,6255 
,4804 
.0232 

0. 

-0. 

0. 

-0. 
0. 

5 
.5531 
,3741 
,0693 
,5943 
,0658 

6 
0.4208 

-0.0964 
0.6654 
0.5324 

-0.2034 

-0. 
-0, 

0, 
-0. 

0. 

7 
,0000 
,0000 
.0000 
.0000 
.0000 

-0, 
-0, 

0, 
-0, 

0, 

P 
,0000 
,0000 
.0000 
,0000 
,0000 

-0, 
-0. 

0. 
-0. 

0. 

9 
,0000 

.oooo 

.0000 

.oooo 

.0000 

LUMPED  PARAMETER  METHOD 
X  COQRDIIN) 


Y  C0QRQ1 IN) 


0.0  120.00000 

360.00000  192.00000 

720.00000  120.00000 

720.00000  0.0 
PINNED  PINNED  END  CONDITION 
ROTATORY  INERTIA  IS  NOT  INCLUDED 

NATUFAL  FREQUENCIESIRAD/SEC) 


NUMBER  OF  SECTIONS 

6.00000 
13.00000 
13.00000 

6.00000 


SPECIFIC  WTILB/CU  IN) 
0.28300 
0.28300 
0.28300 
0.28300 


ARFAISO     IN) 

4.62000 
29.40000 
29.40000 

5.88000 


1 

1 

15, 

.285 

2 
25.365 

3 

75, 

162 

4 
147.574 

5 
269, 

188 

NORMAL IZED  X 

DISPLACEMENTS 

1 

2 
3 

4 
5 

1 
0.0 
0.0 
0.0 
0.0 
0.0 

2 
-0.2987 
0.0689 
-0.2048 
-0. 1023 
-0.7835 

3 
-0.4297 

0.0046 
-0. 1150 

0.  1279 
-0.0052 

-0, 

-0, 

0, 

0, 

-0. 

4 
.4151 
.0936 
.0148 
.22  34 
.0062 

-0, 

-0, 

-0. 

0, 

0. 

5 

3958 
1264 
1238 
0081 
.0015 

6 
-0.3956 
-0. 1783 

0.0013 
-0.2152 

0.0243 

-0. 
-0, 
-0. 
-0, 
0. 

7 
,  '618 
,2574 
.1327 
,1115 
.0081 

-0. 

-0, 

-0, 

0, 

-0. 

p 
,2169 
,2140 
,2013 
.0991 
,6146 

0, 

-0, 

0, 

0, 

-0, 

9 

,0000 
,0000 
,0000 
,0000 
,0000 

tiORMAL 

(2E0  Y 

DISPLACEMENTS 

1 

2 

3 
4 
5 

1 
0.0 
0.0 
0.0 
0.0 
0.0 

2 

0.0 
0.0 
0.0 
0.0 
0.0 

3 
0.0 
0.0 
0.0 
0.0 
0.0 

-0, 

0, 

-0, 

-0. 

0 

4 
,0733 
.4913 
.6488 
.4775 
.0054 

-0. 
0, 
0, 
0, 

-0. 

5 
1698 
6550 

,0444 
5988 

,0332 

6 

-0. 1689 
0.3953 
0.6698 

-0.5178 
0.0807 

0, 
-0, 
-0. 

0. 
-0, 

7 
,0000 
,0000 
,0000 
,0000 
,0000 

0, 
-0. 
-0, 

0, 
-0, 

a 

.0000 
,0000 
,0000 
,0000 
,0000 

0 
-0, 
-0, 

0, 
-0, 

9 
,0000 
,0000 
,0000 
,0000 
,0000 

MOMFNTI  IN     4) 

30.30000 

625.00000 

62^.30000 

°2.  50000 
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LUMPED  PARAMETER  METHOO 
X  COORD! INI 


Y    COQPQI INI 


0.0  120.00000 

360.00000  192.00000 

720.00000  120.00000 

720.00000  0.0 

CLAMPED    CLAMPED    END  CONDITION 

ROTATORY     INERTIA     IS  NOT     INCLUDED 


NUMBER  OF  SECTIONS 

6.00000 
12.00000 
12.00000 

6.00000 


SPECIFIC  WT(LR/ril  INI 
0.28300 
0.28300 
0.28300 
0.28300 


ARFAISO  TN| 

5.08000 
29.40000 
29.40000 

4.62000 


"TMCNTt  IN    41 
f. 50000 
625.00000 

625.00000 
^o.  mnoi 


NATURAL  FREQUENCIEStRAO/SFCI 


NORMALIZED  X  DISPLACEMENTS 


1 


0.0 
0.0 
0.0 
0.0 
0.0 


2 

3 

4 

5 

6 

7 

8 

9 

0, 

.0352 

0, 

,  1452 

0. 

,2293 

0. 

2559 

0. 

,3001 

0, 

.3666 

0. 

,2043 

0 

.0000 

-0. 

.1970 

-0. 

.4195 

-0, 

.4002 

-0. 

,3447 

-0. 

,  3484 

-0, 

.2699 

-0. 

,0935 

-0, 

.0000 

-0 

.1625 

-0, 

,2001 

-0, 

.0672 

-0. 

,1861 

-0, 

,0473 

-0. 

.1721 

-0. 

.1635 

-0, 

.0000 

-0. 

.0609 

0. 

,  1062 

-  0. 

.2120 

-0. 

0137 

-0. 

.2289 

-0, 

.1335 

0. 

,0567 

0 

.0000 

-0. 

.6215 

-0. 

.0252 

0, 

,0147 

-0, 

0385 

-0, 

.0470 

-0, 

.0517 

-0. 

,7495 

-0, 

.0000 

NORMALIZED  Y  DISPLACEMENTS 


1 


0.0 
0.0 

0.0 
0.0 
0.0 


0.0 
0.0 
0.0 
0.0 
0.0 


4 

5 

6 

7 

8 

Q 

-0, 

,4201 

-0. 

.5535 

-0, 

,3325 

-0. 

.0000 

-0 

.0000 

-0 

.0000 

-0. 

,0967 

-0, 

.3741 

-0. 

,3924 

0. 

.0000 

0, 

.0000 

0. 

.0000 

-0. 

,6644 

-0. 

.0700 

0. 

,6241 

-0. 

.0000 

-0, 

.0000 

-0, 

.0000 

-0. 

.5291 

0, 

.5993 

-0. 

,4769 

0, 

.0000 

0 

.0000 

0, 

,0000 

-0. 

,  1999 

0, 

.0662 

0. 

.0234 

-0, 

.0000 

-0, 

.0000 

0. 

.0000 

LUMPEO  PARAMETER  METHOD 
X  CUORDdNI 

0.0  120.00000 

360.00000  192.00000 

720.00000  120.00000 

720.00000  0.0 

CLAMPED  CLAMPEO  END  CONDITION 

ROTATORY  INERTIA  IS  NOT  INCLUDED 

NATURAL  FREOUENCIESIRAD/SEC) 


NUMBER  OF  SECTIONS 

4.00000 

13.00000 

13.00000 

4.00000 


SPECIFIC  WTILB/CU  IN] 
0.28300 
0.28300 
0.28300 
0.28300 


ARFMSO     IN| 

5.88000 
29.40000 
29.40000 

4.62000 


MOMCNT( t«     4) 

82. =0000 

6'5.D0000 

625. OOOOO 

30.30000 


1 

1 
26.180 

35.301 

3 
79. 

662 

4 
148.348 

5 
271, 

,730 

NORMALIZED  X 

DISPLACEMENTS 

1 

1  0.0 

2  0.0 

3  0.0 

4  0.0 

5  0.0 

2 

0.0338 
-0.1956 
-0.1617 

0.0600 
-0.6208 

3 
0. 1427 
-0.4189 
-0. 2011 
-0. 1049 
-0.0240 

0. 
-0. 
-0, 
-0. 

0, 

4 
,2273 
.4009 
.0683 
,2116 
.0174 

0. 
-0. 

-0. 
0. 

-0. 

5 
2538 
3446 
1869 
01  38 
0376 

6 

0.2983 
-0.3494 
-0.0480 

0.2286 
-0.0467 

0, 
-0, 
-0, 

0, 

-0, 

7 
.3650 
.2703 
.1727 
.1326 
.0511 

0, 
-0. 
-0, 

-0, 
-0 

8 
.2032 
.0933 
.1627 
.0558 
,7480 

0 

0, 

0, 

-0, 

-0, 

9 
,0 

.0000 
.0000 
.0000 
.0000 

NORMAL  I ZEO  Y 

DISPLACEMENTS 

1 

1  0.0 

2  0.0 

3  0.0 

4  0.0 

5  0.0 

2 
0.0 
0.0 
0.0 
0.0 
0.0 

3 

0.0 
0.0 
0.0 
0.0 
0.0 

-0, 
-0. 

-0, 

0, 

-0, 

4 
.4230 
,0899 
.6642 
,5333 
.2072 

-0. 
-0. 
-0. 
-0. 
0. 

5 
,5558 
3716 
0709 
,5937 
,0677 

6 

-0.3332 

-0.3959 

0.6237 

0.4802 

0.0219 

-0, 
0, 
-0, 
-0, 
-0. 

7 
,0000 
,0000 
,0000 
,0000 
,0000 

-0, 

0, 

-0, 

-0, 

-0, 

,0000 
,0000 
.0000 
.0000 
,0000 

-0, 
0, 
-0, 
-0, 
-0, 

Q 
.0000 
.0000 
.0000 
.0000 
.0000 
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APPENDIX  C 

FINITE  ELEMENT  PROGRAM 
General 

Program  FINITE  is  a  double  precision  IBM  FORTRAN  IV 

7 
language   digital  computer  program  designed  to  find  the 

principal  vibrational  modes  of  straight  sided  plane  frames 
of  uniform  cross  section.   The  program  is  written  for  the 
IBM  360/67  facility  installed  at  the  Naval  Postgraduate 
School.   The  program  uses  a  f inite-elennent  approach  to 
calculate  the  in-plane  undamped  natural  frequencies  and 
normalized  displacements  of  structures.   The  program  neg- 
lects axial  and  shear  deformation,  but  rotatory  inertia 
may  be  included  if  desired. 

The  program  is  dimensioned  to  handle  a  structure  with 
a  maximum  of  20  members,  and  a  maximum  total  of  43  elements, 
The  program  with  this  dimensioning  requires  total  storage 
of  280K  bytes.   The  properties  of  each  straight  member  must 
be  uniform,  but  may  vary  from  member  to  member. 

The  more  severely  constrained  end  of  the  structure  is 
considered  the  left  end,  and  must  be  pinned  or  fixed.  The 
right  end  may  be  fixed,  pinned,  simply  supported,  or  free. 
The  program  cannot  handle  any  intermediate  supports  or  ex- 
ternally applied  forces. 
Program  Structure 

The  program  consists  of  a  main  body  and  four  subrou- 
tines.  The  program  requires  the  use  of  double  precision 
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sine,  cosine,  arctangent,  absolute  value,  and  square  root 
routines.   If  they  are  not  included  in  the  FORTRAN  library 
they  must  be  provided  by  other  means. 

The  sections  of  the  main  body  and  their  functions  are 
explained  by  comment  cards  in  the  program  listing.   There 
is  also  a  large  leading  section  giving  input  format  infor- 
mation. 

The  subroutine  ELEM  generates  the  elemental  stiffness 
and  mass  matrices  for  use  in  the  main  program. 

The  subroutine  CORNER  performs  the  corner  transforma- 
tion on  the  last  element  to  the  left  of  the  corner  and  the 
first  element  to  the  right  of  the  corner. 

The  subroutines  MAPRIN  and  JACVAT  are  the  same  as 
used  in  the  LUMPED  program.   These  subroutines  perform  the 
same  function  in  the  FINITE  program  as  they  did  in  the 
LUMPED  program.   Once  again  the  rows  and  columns  of  the 
mass  and  stiffness  matrices  deleted  in  applying  end  condi- 
tions are  replaced  by  zeroes.   In  the  mass  matrix,  each 
diagonal  element  corresponding  to  a  deleted  row  and  column 
is  replaced  with  unity. 
Input 

The  input  required  by  the  FINITE  program  is  the  same 
as  that  needed  for  the  LUMPED  program  with  one  exception. 
The  number  of  elements  into  which  a  member  is  subdivided 
must  be  an  even  number.   The  input  information  for  the 
LUMPED  program  is  found  on  page  5  8  of  Appendix  B. 
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Output 

The  form  and  type  of  output  for  the  FINITE  program 
is  identical  to  that  of  the  LUMPED  program  and  is  discussed 
on  page  59  of  Appendix  B.   In  order  to  specify  the  number 
of  frequencies  desired,  card  number  999  in  the  program  list- 
ing is  used.   The  number  of  frequencies  may  not  exceed  twice 
the  number  of  elements  used. 
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FINITE  ELEMENT  METHOD 

X  COORDIINI  Y  CriORQ(IN)         NUMBER  OF  SECTIONS   SPECIFIC  WTILR/CU  INI     ARFAISO  INI        MOMFNTIIN  41 

2.00000  0.28300  4.62000  30.30000 

0.0  120.00000 

4.00000  0.28300  29.40000  675.00000 

360.00000  192.00000 

4.00000  0.28300  20.40000  6»5.D000O 

720.00000  120.00000 

2.00000  0.28300  5.88000  82.50000 

720.00000  0.0 

CLAMPED  CLAMPED  END  CONDITION 

ROTATORY  INERTIA  IS  NOT  INCLUDED 

NATURAL  FREQUENCIESIRAO/SECI 


1 

1 
25.752 

2 
34.679 

3 
79. 

,647 

4 
150.444 

5 
280. 

,597 

NORMALIZED  X 

DISPLACEMENTS 

1 

1  0.0 

2  0.0 

3  0.0 

4  0.0 

5  0.0 

2 
0.2051 
0.0932 
0.1615 
0.0551 
-0.7539 

3 

0. 3680 

0.2685 

0. 1677 

-0. 1330 

-0.0475 

0. 

0. 

0. 
-0. 
-0, 

4 
.3016 
,3468 
,0421 
.2284 
.0431 

0. 

0. 

0, 
-0. 
-0. 

5 
.2578 
,  3440 
.1812 
,0126 
,0348 

6 
0.2313 
0.3988 
0.0612 
0.2130 
0.0153 

0, 
0. 
0, 
0. 
-0, 

7 
.1476 
.4196 
.1946 
.1077 
,0222 

0, 

0. 

0. 

-0, 

-0, 

8 
.0365 
.  1980 
.1604 
.0587 
.6213 

0, 
0, 
0, 
0, 
0, 

,  0 
.0 
,  0 
,0 
,0 

NORMAL  1  ZED  Y 

DISPLACEMENTS 

1 

1  0.0 

2  0.0 

3  0.0 
*      0.0 
5      0.0 

2 

0.0 
0.0 
0.0 
0.0 
0.0 

3 
0.0 
0.0 
0.0 
0.0 
0.0 

0. 
-0, 

0. 

0, 
-0. 

4 
,3317 
,3912 
,62  79 
,4771 
,0217 

0. 
-0. 
-0. 
-0. 
-0, 

5 
,5509 
,3775 
,0674 
,6016 
.0632 

6 

0.4183 

-0.1036 

-0.6671 

0.5263 

0.1872 

0. 
0, 
0, 
0. 
-0, 

7 
.0000 
,0000 
.0000 
,0000 
,0000 

0, 
0, 
0. 
0. 
0, 

8 
.0 
,0 
.0 
,0 
.0 

0, 
0. 
0, 
0, 

0. 

.0 
,0 
,0 

,0 
,0 

FINITE  ELEMENT  METHOD 

X  COORDIIN)           Y  COQROIINI         NUMBER  OF  SECTIONS  SPECIFIC  WTILB/CU  IN)      AREAISQ  INI  MnMFNT<TW  41 

2.00000              0.28300              4.62000  3O.30O00 

4.00000              0.28300             29.40000  675.00000 

4.00000              0.28300             29.40000  675.00010 

2.00000               0.28300               5.88000  «2.5000O 


0.0  120.00000 

360.00000  192.00000 


720.00000  120.00000 

720.00000  0.0 


PINNED  PINNEO  END  CONDITION 
ROTATORY  INERTIA  IS  NOT  INCLUDED 

NATURAL  FREQUENCIESIRAO/SECI 


1 

1 
15. 

256 

2 
25.364 

3 
75, 

,667 

4 
149.820 

5 
278. 

140 

NORMALIZED  X 

DISPLACEMENTS 

ooooo 
ooooo 

2 

0.2982 

0.0689 

-0.2040 

-0.0°93 

-0.7862 

3 

0.4299 
0.0054 

-0. 1144 
0. 1789 

-0.0051 

4 
0.4152 

-0.0927 
0.0154 
0.2238 

-0.0060 

0, 
-0. 
-0, 

0. 
0. 

5 
,3960 
,1257 
,1232 

,0078 
,0013 

6 

0.3957 
-0.1777 

0.0021 
-0.2159 

0.0229 

0, 
-0, 
-0, 
-0, 

0, 

7 
.3672 
,2569 
,1319 
,1133 
.0078 

0. 
-0 
-0, 

0, 
-0, 

0 
.2170 
,2132 
.1993 
.0061 
,6119 

0, 
0. 
0. 

0, 
0, 

,0 
,0 
,0 
.0 
,0 

NORMALIZED  Y 

DISPLACEMENTS 

1 

1  0.0 

2  0.0 

3  0.0 

4  0.0 

5  0.0 

2 
0.0 
0.0 
0.0 
0.0 
0.0 

3 
0.0 
0.0 
0.0 
0.0 
0.0 

4 

0.0733 

0.4906 

-0.6493 

-0.4743 

0.0046 

0. 
0. 
0. 
0. 
-0, 

5 
,  1693 
,6559 
,0438 
,6054 
.0372 

6 

0. 1678 
0.3961 
0.6702 
-0.5130 
0.0756 

0, 
-0, 
-0, 
-0, 

0, 

7 
,0000 
.0000 
.0000 
.0000 
.0000 

0, 
0 
0, 
0 
0, 

8 
.0 
.0 
.0 
.0 
,0 

0, 
0, 
0, 
0, 
0, 

.0 
,0 

.0 
.0 
,0 
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FINITE  ELEMENT  METHOD 

X    COORO(IN)  Y    CnpRDIIN)  NUMBER    OF    SECTIONS       SPFCIFIC    WTILR/CII    IN)  ARFAISO     INI  MflMFNTMN    4) 


0.0  120.00000 

360.00000  192.00000 

720.00000  120.00000 

720.00000  0.0 

CLAMPEO    CLAMPED    END  CONDITION 

ROTATORY    INERTIA    IS  NOT    INCLUDED 

NATURAL    FREQUENCIESIRAO/SEC) 


2.00000  0.28300  5.PR000  82.*Onoo 

2.00000  0.28300  29.40000  6?"!. OOOOO 

2.00000  0.28300  29.40000  62*. OOOOO 

2.00000  0.28100  4.62000  30.30000 


1 

1 
25.761 

2 

34.694 

3 
79, 

,974 

4 
152.647 

5 
311, 

542 

NORMALIZED  X 

DISPLACEMENTS 

1 

1  0.0 

2  0.0 

3  0.0 

4  0.0 

5  0.0 

2 
-0.0366 

0.1980 
-0.1596 

0.0576 
-0.6113 

3 
-0.1478 
0.4194 
-0.1928 
-0. 1079 
-0.0122 

-0, 

0, 

-0, 

-0. 

0. 

4 
.2315 
.3985 
,0592 
,2134 
,0212 

5 

-0.2580 
0.3437 

-0.1793 
0.0121 

-0.0237 

6 
-0.3018 

0.3463 
-0.0401 

0.2278 
-0.0320 

-0, 
0, 

-0, 
0, 

-0, 

7 
.3681 
,2680 
,  1659 
,1321 
.0351 

-0, 
0, 
-0, 
-0, 
-0, 

8 
,2051 
,0979 
.1608 
,0545 
,7690 

0 
0 
0 
0 
0, 

9 

.0 

,n 

.0 
,0 
.0 

NORMALIZED  Y 

DISPLACEMENTS 

1  0.0 

2  0.0 

3  0.0 

4  0.0 

5  0.0 

2 

0.0 
0.0 
0.0 
0.0 
0.0 

3 
0.0 
0.0 
0.0 
0.0 
0.0 

0, 
0. 

-0. 
0. 

-0. 

4 
.4182 
.1044 
.66  82 
.5272 
.1670 

5 

0.5507 

0.3785 

-0.0673 

-0.6001 

0.0573 

6 
0.3315 
0.3916 
0.6291 
0.4786 
0.0157 

-0, 
0. 

-0, 

0, 

-0. 

7 
.0000 
,0000 
,0000 
.0000 
,0000 

0, 
0, 
0, 
0. 
0, 

8 
,0 
,0 
,0 
.0 
,0 

0, 
0, 
0, 
0, 
0, 

9 

.0 
,0 
.0 

.0 
,0 

FINITE    ELEMENT    METHOD 

X    COORD(IN)                           Y    COORD(IN)                      NUMBER    OF    SECTIONS  SPECIFIC    WTIL8/CU    IN)  AOEAISO     I  »l  I  MDMFNTIIN    4) 

4.00000                                   0.28300                                    5.88000  82.50000 

4.00000                                   0.28300  29.40000  625.00000 

4.00000                                      0.28300  29.40000  625.30000 

4.00000                                   0.28300  4.62000  30.30000 


0.0  120.00000 

360.00000  192.00000 


720.00000  120.00000 

720.00000  0.0 


CLA«PED    CLAMPED    ENO    CONOITION 
ROTATORY     INERTIA     IS    NOT     INCLIIOED 

NATURAL  FREOUENCIESIRAD/SEC) 


1 

1 
25.752 

2 
34.679 

3 
79, 

,647 

4 
150.444 

5 
280. 

,594 

►WKMALIZED  X 

DISPLACEMENTS 

1 

1  0.0 

2  0.0 

3  0.0 

4  0.0 

5  0.0 

2 

-0.0365 
0. 1980 

-0. 1604 
0.0587 
0.6211 

3 
-0. 1476 

0.4196 
-0. 1946 
-0.1077 

0.0221 

4 
-0.2313 

0.3988 
-0.0612 
-0.2130 
-0.0153 

-0, 
0. 

-0, 
0. 

0, 

5 
2578 
3440 
,1812 
0126 
0347 

6 

-0.3016 
0.3468 

-0.0421 
0.2284 
0.0430 

-0. 
0. 

-0, 
0, 
0, 

7 
.3680 
.2685 
.1677 
.1330 
.0474 

-0, 

0, 

-0, 

-0, 

0, 

a 
,2051 
.091? 
.1615 
.055? 
,7541 

0, 
0. 
0, 
0, 
0, 

.0 
,0 
,0 
.0 
.0 

NORMALIZED  Y 

DISPLACEMENTS 

1 

1  0.0 

2  0.0 

3  0.0 

4  0.0 

5  0.0 

2 

0.0 
0.0 
0.0 
0.0 
0.0 

3 
0.0 
0.0 
0.0 
0.0 
0.0 

4 
0.4183 
0.1036 
-0.6671 
0.5263 
0.1871 

0, 
0. 

-0. 

-0. 

-0, 

5 
5509 
3775 
0674 
,6016 
,06  32 

6 
0.3317 
0.3912 
0.6279 
0.4771 
-0.0217 

-0, 
0, 

-0. 
0, 
0. 

7 
.0000 
.0000 
.0000 
.0000 
,0000 

0, 
0. 
0, 
0, 
0, 

8 
,0 
,0 
.0 
,0 
,0 

0, 
0, 
0, 
0, 
0, 

,0 

,0 
,0 

,0 
.0 
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APPENDIX  D 

TRANSFER  MATRIX  PROGRAM 
General 

Program  DYNRES  is  a  double  precision  IBM  FORTRAN  IV 

7 
language   digital  computer  program  designed  to  obtain 

"exact"  solutions  for  the  principal  modes  of  straight 

sided  plane  frames  of  uniform  cross  section.   These  modes 

are  used  as  a  standard  for  comparing  the  LUMPED  and  FINITE 

programs.   Axial  and  shear -deformation ,  as  well  as  rotatory 

inertia,  are  neglected. 

9 

The  program  uses  the  method  of  false  position   to 

assist  convergence  to  the  frequencies.   Convergence  is  ob- 
tained when  either  the  absolute  value  of  the  frequency  de- 

—  16 
terminant  is  less  than  10    ,  or  the  frequency  is  bracketed 

-9 

such  that  the  interval  is  less  than  10 

The  program  is  dimensioned  for  use  with  structures 
having  no  more  than  nine  members.   The  program  requires  a 
total  storage  of  50K  bytes.   Each  member  must  have  uniform 
properties ,  but  the  properties  may  vary  from  member  to  mem- 
ber. 

The  more  severely  constrained  end  of  the  structure  is 
designated  the  left  end.   This  end  must  be  either  fixed  or 
pinned.   The  right  end  is  allowed  any  of  the  four  common 
boundary  conditions:   fixed,  pinned,  simply  supported,  or 
free.   The  structure  may  not  have  any  intermediate   supports 
or  externally  applied  forces. 
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Program  Structure 

The  program  consists  of  the  main  body  and  two  subrou- 
tines.  The  program  requires  the  use  of  double  precision 
sine,  cosine,  arctangent,  hyperbolic  sine,  hyperbolic  co- 
sine, and  square  root  routines.   If  these  routines  are  not 
available  in  the  FORTRAN  library,  they  must  be  supplied  by 
other  means.   The  main  body  is  prefaced  with  a  comment  sec- 
tion containing  the'  format  and  makeup  of  the  input  deck. 
The  functions  of  the  various  sections  of  the  main  program 
are  explained  in  comment  cards  spaced  throughout  it. 

The  subroutine  EVAL  generates  the  transfer  matrices 
and  evaluates  the  frequency  determinant  after  a  trial  fre- 
quency has  been  determined  by  the  main  program. 

The  subroutine  MAPRIN  is  the  output  subroutine.   It 
performs  in  the  same  manner  as  discussed  in  the  LUMPED 
program  on  page  5  8  of  Appendix  B. 
Input 

The  input  required  by  the  DYNRES  program  is  basically 
the  same  as  for  the  LUMPED  program  (page  5  8  of  Appendix  B) . 
The  DYNRES  program  requires  the  coordinates  of  the  midpoint 
of  each  member  as  well  as  the  end  and  corner  coordinates. 
No  further  subdivision  of  the  frame  is  required.   The  rest 
of  the  data  is  the  same.   Format  and  order  of  data  cards  are 
given  at  the  front  of  the  main  program. 
Output 

The  coordinates  of  ends,  midpoints,  and  corners  are 
automatically  printed  by  the  program.   The  properties 
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(modulus  of  elasticity,  specific  weight,  area,  moment  of 
inertia)  between  each  set  of  coordinates  are  also  printed 
The  rest  of  the  output  is  identical  with  that  of  the 
LUMPED  program.   The  number  of  frequencies  solved  for  is 
specified  by  a  control  card.   Its  makeup  is  shown  in  the 
comment  section  preceding  the  program. 
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TRANSFER  MATRIX  METHOD 

X  COORD*  IN)  Y  COORD!  IN) 

0.0  60.00000 

0.0  120.00000 

180.00000  156.00000 

360.00000  192.00000 

540.00000  156.00000 

720. OOOOO  120.00000 

720.00000  60.00000 

720.00000  0.0 
CLAMPEO  CLAMPEO  END  CONDITION 

ROTATORY  INERTIA  IS  NOT  INCLUDED 


MODULUS  OF  ELASTICITYIPSI) 
30000000.00000 
300*0000. OOOOO 
30000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 


SPECIFIC  WT(LB/CII  IN)   AREAISO  IN) 

MOMFNTI IN  4) 

0.28300 

4.62000 

■»0.3000" 

0.28300 

4.62000 

30.30000 

0.28300 

29.40000 

625.00000 

0.28300 

29.40000 

625.00000 

0.28300 

29. 40000 

625.00000 

0.28300 

29.40000 

625.00000 

0.28300 

5.88000 

82.50000 

0.28300 

5.08000 

02.50000 

NATURAL  FREOUENCIESIRAO/SEC) 


NORMALIZED  X  DISPLACEMENTS 


0.0 
0.0 
0.0 
0.0 
0.0 


2 
0.2051 
-0.0932 
0.1616 
0.0551 
0.7538 


3 

0.3679 
-0.2686 

0.  1678 
-0.1330 

0.0479 


0.3015 
-0.3468 

0.0422 
-0.2284 

0.0435 


5 

0.2577 
-0.3441 

0.1813 
-0.0127 

0.0352 


6 
2311 
3988 

0614 
2129 

0149 


7 
0.1475 
-0.4196 
0.1948 
0.1077 
0.0225 


8 

0.0'".4 
-0.1981 

0.1604 
-0.0587 

0.6213 


-0.0000 
-0.0000 
-0.0000 
-0.0000 
0.0000 


NORMALIZED  Y  DISPLACEMENTS 


1 

0.0 

2 

0.0 

3 

0.0 

4 

0.0 

5 

0.0 

-0.0000 

o.oooo 

-0.0000 
-0.0000 
-0.0000 


-0.0000 
0.0000 

-0.0000 
0.0000 

-0.0000 


3318 
3911 
62  79 
4770 


5 

0.5511 

0.3775 

-0.0674 

-0.6017 

0.0633 


0.4183 
0.1036 

-0.6670 
0.5262 

-0.1874 


-0.0000 
-0.0000 
-0.0000 
-0.0000 
0.0000 


-0.0000 
-0.0000 
-0.0000 
-0.0000 
0.0000 


-0.0000 
-0.0000 
-0.0000 
-0.0000 
0.0000 


TRANSFER  MATRIX  METHOD 

X  COORD! IN)  Y  COORC(IN) 

0.0  60.00000 

0.0  120.00000 

180.00000  156.00000 

360.00000  192.00000 

540.00000  156.00000 

720.00000  120.00000 

720.00000  60.00000 

720.00000  ..  0.0 
PINNED  PINNEO  END  CONDITION 

ROTATORY  INERTIA  IS  NOT  INCLUDED 

NATURAL  FREOUENCIESIRAO/SEC) 


"OQULUS  OF  ELASTICITYIPSI) 
30000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 
30000000.00000 


SPFCI FIC  UTIL B/CU  INI 

AREAISO  IN) 

MOMENT!  IN  4) 

0.28300 

4.62000 

30.30000 

0.28300 

4.6'000 

30.30000 

0.28300 

29.40000 

625. OOOOO 

0.28300 

29.40000 

6'5. OOOOO 

0.28300 

29.40000 

625.00000 

0.28300 

2Q.40000 

625.00000 

0.28300 

5.88000 

02.50000 

0.28300 

5.88000 

82. ^0000 

1 

15, 

.257 

25.365 

•   75, 

.652 

149.662 

277 

.069 

NORMALIZED  X 

DISPLACEMENTS 

1 

1  0.0 

2  0.0 

3  0.0 

4  0.0 

5  0.0. 

2 
0.2982 
0.0689 
0.2040 
0.0993 
0.7861 

3 
0.4299 
0.0064 
0.1145 
-0.1289 
0.0052 

4 

0.4152 

" -0.0927 

-0.0153 

-0.2238 

0.0061 

0. 
-0. 

0. 
-0. 

-0, 

5 
,3960 
,1257 
.1233 

.0078 
,0013 

6 

0.3957 
-0.1777 
-0.0020 

0.2159 
-0.0229 

0, 

-0. 

0. 

0, 
-0. 

7 
.3622 
.2569 
.1320 
,1133 
,0077 

0. 
-0, 

0, 
-0, 

0, 

8 
.2170 
,2132 
,1993 
.0961 
.6120 

-0, 
0 
-0, 
-0 
-0, 

9 

.0000 
,0000 
,0000 
,0000 
,0001 

NORMALIZED  Y 

DISPLACEMENTS 

1 

1  0.0 

2  0.0 

3  0.0 

4  0.0 

5  0.0 

2 

-0.0000 
-0.0000 
-0.0000 
-0.0000 
-0.0000 

3 
-0.0000 
-0.0000 
-0.0000 
0.0000 
-0.0000 

•   4 

. 0.0733 

*   0.4906 

•   0'.6493 

0.4743 

-0.0046 

0. 

0, 
-0, 
-0. 

0. 

5 
,1693 
,6559 
,0438 
,6055 
.0323 

6 
0.1678 
0.3961 

-0.6702 
0.5130 

-0.0758 

-0, 
0. 
-0. 
-0. 
-0. 

7 
,0000 
,0000 
,0000 
,0000 
.0000 

-0. 
0. 
-0. 
-0, 
-0. 

8 
,0000 
,0000 
,0000 
,0000 
,0000 

-0, 

0, 
-0, 
-0, 

-0, 

9 

,0000 
,0000 
,0000 
,0000 
,0000 
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